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ABSTRACT  OF  THE  DISSERTATION 


Nonlinear  Acoustics  in  a  Dispersive  Continuum: 

Random  Waves,  Radiation  Pressure,  and  Quantum  Noise 

by 

Michael  Anthony  Cabot 
Doctor  of  Philosophy  in  Physics 
University  of  California,  Los  Angeles,  1983 
Professor  Seth  J.  Putterman,  Chair 

The  nonlinear  interaction  of  sound  with  sound  is  studied  using 
dispersive  hydrodynamics  which  is  derived  from  a  variational  princi¬ 
ple  and  the  assumption  that  the  internal  energy  density  depends  on 
gradients  of  the  mass  density.  The  attenuation  of  sound  due  to  non¬ 
linear  interaction  with  a  background  is  calculated  and  is  shown  to 
be  sensitive  to  both  the  nature  of  the  dispersion  and  decay  band- 
widths.  The  theoretical  results  are  compared  to  those  of  low  tempera¬ 
ture  helium  experiments.  A  kinetic  equation  which  describes  the  non¬ 
linear  self-interaction  of  a  background  is  derived.  When  a  Debye-type 
cutoff  is  imposed,  a  white  noise  distribution  is  shown  to  be  a  sta¬ 
tionary  distribution  of  the  kinetic  equation.  Zero  point  motion  is 
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introduced  into  the  classical  hydrodynamics  through  a  renormalization 
scheme,  which  imposes  the  requirement  that  in  a  sound  scattering  the 
zero  point  motion  does  not  lose  energy.  The  form  of  the  zero  point 
motion  is  then  determined  by  the  kinetic  equation  and  a  derived  fluid 
law  that  is  analogous  to  Wien's  displacement  law  for  electromagnetic 
radiation.  The  kinetic  equation  with  zero  point  motion  included  is 
shown  to  have  a  Planck  distribution  as  a  stationary  solution.  An 
H-theorem  is  presented.  The  attenuation  and  spectrum  of  decay  of  a 
sound  wave  due  to  nonlinear  interaction  with  zero  point  motion  is 
calculated.  In  one  dimension,  the  dispersive  hydrodynamic  equations 
are  used  to  calculate  the  Langevin  and  Rayleigh  radiation  pressures 
of  wave  packets  and  solitary  waves. 


INTRODUCTION 


The  hydrodynamical  equations  for  a  classical  one  component  fluid 
constitute  a  contracted  description  in  which  knowledge  of  a  complete 
set  of  five  field  variables  (e.g.  p,s,v;  the  density,  specific  entropy, 
and  velocity)  uniquely  specifies  the  state. ^  These  equations  are  in 
general  nonlinear  due  to  the  convective  terms,  which  arise  as  a  con¬ 
sequence  of  the  requirements  imposed  by  Qalllean  relativity,  as  well  as 
the  form  of  the  equation  of  state  for  pressure  "?"  as  a  function  of 
density  which  is  needed  to  close  the  theory.  The  equations  of  hydro¬ 
dynamics  are  a  unique  consequence  of  the  assumption  of  a  local  descrip¬ 
tion  in  terms  of  five  field  variables  as  a  complete  set,  and  Pascal’s 
law  for  the  isotropy  of  the  stresses  (i.e.  the  stress  tensor  in  the 
rest  frame  must  be  diagonal) . 

For  disturbances  of  sufficiently  small  amplitude  the  hydrodynamics 
can  be  linearized  in  the  amplitude  of  deviation  from  the  equilibrium 
state.  One  finds  that  there  exist  traveling  wave  solutions  (sound 

waves)  where  the  frequency  w  is  related  to  the  wave  number  k  by 

m  L 

u  ■  cit  where  c  »  (—  ]2  . 

Such  a  linear  relation  is  said  to  imply  no  dispersion.  On  the  linear 
level  superposition  applies,  sound  waves  run  through  each  other  without 
interacting  so  that  the  propagation  of  one  wave  is  unaffected  by  the 
presence  of  other  waves.  Nonlinearities  cause  a  propagating  wave  to 
interact  with  other  waves  so  as  to  scatter  energy  out  of  the  original 
channel  and  with  Itself  so  as  to  distort  and  eventually  form  a  shock- 
front.  The  main  goal  of  this  thesis  is  to  discuss  the  consequences 


for  sound  propagation  of  including  anisotropic  stresses  (and  therefore 
dispersion)  as  well  as  nonlinear  effects. 


The  nonlinear  dispersive  equations  are  derived  from  a  variational 

2 

principle  (for  an  energy  density  which  is  a  function  of  (Vp)  )  in 
Chapter  I,  where  their  general  properties  (e.g.  Kelvin,  Bernoulli  and 
conservation  laws)  are  discussed.  At  lowest  order  the  nonlinear 
equations  take  the  form: 

2  2 

Op  «  GV  p  (0-1) 


where 


2  2  4 

c  7  +  a7 


and  G  is  a  nonlinear  coupling  coefficient.  This  equation  will  always 
be  solved  in  a  perturbation  approach.  For  instance  considering  a  super¬ 
position  of  two  waves  which  satisfy  the  lowest  order  equation  []p*0; 
so  that 

p  ■  picos(ki*r  -  ui^t)  ♦  P2Cos(k2'r  -  u>2t)  . 


one  finds  that  the  leading  order  correction  to  0  due  to  the  interaction 
of  the  two  waves  is  p ^  where 

0pi2  ■  GpiP2^  l  cos[(ki  t  k2)'r  -  (uj  ±  u»2)t] 

± 


so  that  nonlinearities 


create  sum  and  difference  frequencies. 


2 
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In  Appendix  A  we  consider  Che  energy  and  momentum  of  a  ID  wave 
packet  localized  in  wave  number  and  propagating  according  to  (0-1) .  As 
a  result  of  the  nonlinearities  some  of  the  energy  and  momentum  are 
contained  at  (almost)  zero  wavelength  and  therefore  part  of  the  energy 

and  momentum  propagate  at  a  different  velocity  than  the  original  packet 

3 

if  the  fluid  is  dispersive.  The  various  radiation  pressures  of  these 
packets  are  calculated  and  compared  with  the  properties  of  solitons. 

As  a  propagating  wave  interacts  with  a  background  of  sound  the 
waves  created  at  the  sum  and  difference  frequencies  takes  energy  from 
the  imposed  wave  and  accounts  for  an  attenuation  which  was  calculated 
by  Westervelt^  for  the  nondispersive  cane  (ct*0) .  The  main  contribution 
to  the  scattering  (attenuation)  arises  from  terms  on  the  right  hand 
side  of  (0-1)  which  are  at  resonance  so  that 

-  2  j. 

<*»1  ±  o»2  -  c  ±  k2|  +  a  Iki  ±  k2  I  . 

In  the  limit  of  no  dispersion  this  occurs  when  two  waves  are  parallel 

-►  ■  i  Hh 

(k^( Ikj).  The  presence  of  a  small  reversible  dispersion  can  have  a 
dramatic  effect  since  for  a>0  there  is  a  resonance  possible  at  some 
finite  angle  whereas  for  a<0  no  resonance  is  possible  and  the  scattering 
of  sound  is  dramatically  reduced.  In  Chapter  II  the  scattering  of 
sound  by  sound  for  the  dispersive  medium  is  calculated.  We  find  that 
Westervelt's  formula  is  multiplied  by  a  factor  of  2  for  :x>0  and  by  "0" 
for  a<0.  This  effect  applies,  however,  only  when  the  dispersion  is 
larger  than  the  bandwidth  which  results  from  the  scattering.  In 
Appendix  B  the  detailed  transition  formulas  that  apply  in  this  inter¬ 
mediate  regime  are  calculated  as  well  as  an  extra  (usually  logarithmic) 
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dispersion  of  sound  due  to  scattering  by  sound. 

The  results,  and  especially  the  approach  employed  should  be  of 
relevance  to  dielectric  media  (e.g.  Helium,  Quartz)  at  low  temperature. 
A  dielectric  is  a  medium  in  which  the  electrons  are  not  free  to  flow  so 
that  the  only  motion  possible  is  a  vibration  of  the  constituent  ions  or 
molecules.  These  oscillations  are  of  course  sound  waves  and  in  ad¬ 
dition  to  the  impressed  wave  there  will  be  a  background  of  random 
oscillations  due  to  the  internal  energy  (or  temperature)  of  the  dielec¬ 
tric  which  must  manifest  itself  in  the  form  of  sound  (noise) .  Based 
upon  our  presumption  that  the  nonlinear  continuum  mechanics  can  be 
applied  to  the  sound  wave  as  well  as  the  thermal  noise  we  can  use  these 
results  of  the  macroscopic  continuum  to  calculate  the  attenuation  coef¬ 
ficient  of  sound  in  dielectrics.  Liquid  helium  is  a  particularly  rich 
dielectric  since  a>0  for  P<18  atmospheres  and  less  than  zero  at  higher 
pressure.  In  Chapter  II  some  of  the  results  are  compared  with  the 
measured  ultrasonic  attenuation  Droperties  in  liquid  helium.^ 

While  nonlinear  effects  can  cause  energy  to  be  scattered  out  of 
a  given  channel  they  can  also  lead  to  processes  in  the  background 
whereby  two  waves  collide  and  throw  energy  back  into  the  given  channel. 
These  effects  are  much  like  restituting  collisions  in  kinetic  theory 
and  are  included  in  Chapter  III  where  a  Boltzmann  type  equation^  for 
sound  waves  i3  derived  from  the  hydrodynamics.  The  stationary  solution 
to  the  nonlinear  Boltzmann  equation  for  sound  is  white  noise  (equipar- 
tition)  thus  deriving  the  noise  spectrum  which-  is  normally  imposed  upon 
the  linearized  theory  via  the  fluctuation  dissipation  theorem. 
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All  of  the  results  discussed  above  (both  theory  and  experiment) 
apply  to  the  regime  where  huj/kgT<l.  In  the  strong  quantum  limit 
hui/k  T>1  sound  waves  should  also  attenuate  due  to  a  spontaneous  decay 

D 

g 

as  first  calculated  by  Slonimskii  using  a  quantized  theory.  In 

Chapter  III  we  show  how  these  results  can  also  be  obtained  from  the 

classical  nonlinear  hydrodynamics  by  including  the  possibility  of 

scattering  from  a  zero-temperature  noise  (or  so-called  zero  point 
9 

motion).  By  the  second  law  of  thermodynamics  the  zero-temperature 
state  cannot  lose  energy.  Thus  the  scattering  calculation  must  be 
renormalized  so  as  not  to  remove  energy  from  the  zero -temperature  noise. 
By  selecting  this  configuration  out  of  the  classical  hydrodynamics  we 
generalize  the  Westervelt  equation  to  include  attenuation  due  to  scat¬ 
tering  by  quantum  as  well  as  thermal  noise.  The  Boltzmann  equation  for 
sound  is  also  generalized  so  as  to  include  a  zero  point  noise  defined 
as  above.  The  resulting  equation  for  the  quantum  theory  of  interacting 
sound  waves  has  a  stationary  solution  given  by  Planck’s  law  plus 
harmonic  zero  point  motion.  Planck's  constant  enters  as  an  undetermined 
constant  which  characterizes  the  most  general  solution  to  the  balance 
equation.  The  uniqueness  of  this  solution  as  well  as  the  corresponding 
"H"  theorem  are  discussed. 

A  central  theme  of  this  work  is  that  many  results  normally  thought 
to  be  obtainable  only  through  microscopic  and  or  quantum  calculations 
can  be  interpreted  and  obtained  directly  from  the  classical  nonlinear 
hydrodynamics  plus  the  second  law.  Thus  the  approach  developed  here 
might  be  of  value  in  understanding  quantum  effects  in  complex  N  body 
systems.  In  particular  when  it  is  desired  to  understand  quantum 
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noise  in  macroscopic  systems  such  as  Josephson  junctions^  it  should  be 
realized  that  the  usual  prescription^  for  writing  down  a  quantum 
Langevin  equation  does  not  take  into  account  the  renormalization 
developed  in  Chapter  III. 

As  the  energy  density  of  the  sound  field  obeys  a  Boltzmann  type 

equation  one  can  ask  if  there  are  compressional  modes  of  the  spectral 

intensity  of  sound.  This,  so-called  second  sound  (thermal  waves)  exists 

in  liquid  helium  and  can  be  obtained  via  the  procedure  developed  herein. 

In  this  vein  it  has  also  recently  been  shown  that  the  complete  two 

fluid  hydrodynamics  can  be  deduced  from  the  nonlinear  continuum 

1 2 

mechanics  of  the  one-component  classical  fluid. 
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CHAPTER  I. 


THE  NONLINEAR  DISPERSIVE  CONTINUUM  MECHANICS 

AC  low  temperatures  sound  in  liquid  helium  and  other  dielectrics 
decays  due  to  scattering  from  other  sound  waves  in  the  substance.  We 
regard  these  ocher  waves  as  the  thermal  and  quantum  noise  present  in 
the  material.  As  two  sound  waves  can  interact  only  as  a  result  of  non¬ 
linear  processes  a  description  of  these  effects  must  be  based  upon  the 
complete  nonlinear  hydrodynamics.  Furthermore  the  interaction  of  sound 
waves  is  strongest  when  the  frequencies  and  wave  numbers  are  in 
resonance  or 

<*>(k)  ■  ♦  sij  and  k  •  k;  t  k2  ,  (1-1) 


where  oi(k)  is  the  dispersion  law  for  the  medium,  and  we  have  in  mind  an 
interaction  of  two  waves  with  wavenumbers  k^  and  k2>  Whether  or  not 
the  condition  (1-1)  can  be  satisfied  depends  crucially  upon  the  dis¬ 
persion  in  the  medium  i.e.  on  Che  second  derivative  of  u  with  respect  to 
k  which  divides  into  3  cases: 


2 

d  a)  ,  . 

2  '  0  upward  or  anomalous  dispersion, 
dk 

2 

<J  <j j 

— r  <  0  downward  or  normal  dispersion,  and 
dk 

2 

d  w 

— ^  *  0  nondispersive  (or  so-called  semi-dispersion), 
dk 


(1-2A) 


( 1-2B) 


(1-20 
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For  cases  (1-2A)  and  (1-2C)  resonances  exist  whereas  for  case  (1-2B) 
there  is  no  resonant  configuration  for  the  three  wave  processes  as  de¬ 
scribed  by  Eq.(l-l).  Thus  the  nonlinear  effects  and  therefore  the 
scattering  of  sound  by  sound  is  sensitive  to  dispersion  so  that  a  basic 
description  of  sound  interacting  with  sound  must  be  dispersive  as  well 
as  nonlinear  .  We  will  see  that  for  case  (1-2B)  the  attenuation  is 


zero  to  this  order  whereas  for  case  (1-2A)  the  attenuation  is  twice 

that  given  by  (1-2C)  though  there  is  a  transition  region.  These  three 

cases  will  be  used  as  a  guide  for  understanding  the  attenuation  of  sound 
4 

in  He  at  low  temperature. 

In  the  continuum  mechanical  approach  dispersion  arises  from  an 
internal  energy  which  depends  upon  gradients  of  the  mass  density.  The 
simplest  model  in  which  this  occurs  is  that  of  particles  on  springs  for 
which  the  Lagrangian  is 


1*2  2 
L  ■  2  ?  *  *<TVl”V  ] 


(1-3) 


where  is  the  displacement  of  the  particle  of  mass  m  from  its 
equilibrium  position  "ja"  where  "a”  Is  the  equilibrium  interparticle 
spacing  and  $  is  the  spring  constant.  We  procede  to  the  continuum  limit 


by  setting 


(1-4) 


8 


30  ehac 
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We  see  that  for  this  simple  model  the  potential  energy  terms  which 

depend  upon  gradients  of  the  density  lead  to  the  dispersion. 

To  obtain  the  general  nonlineaT  equations  including  dispersion  we 

apply  Hamilton’s  principle  to  a  fluid  Lagrangian  density  given  by  the 

1  2 

difference  of  kinetic  (y  pv  )  and  internal  energy  density  U  where 


0  -  U(p , s ,Vp) 


(1-10) 


so  that 


dj(ypv2  -  U)drdt  *  0 


(1-11) 


This,  however,  must  be  subject  to  the  constraints  of  mass  conservation 
and  entropy  conservation  which  are  not  contained  in  the  symmetry  of  the 
original  Lagrangian : 


bp  ♦  ♦ 

+  7*pv  «  0  and 
bt 


(1-12) 


+  7-psv  »  0 


(1-13) 


Incorporating  (1-12)  and  (1-13)  in  (1-11)  with  undetermined  Lagrange 
multipliers  $,  8  leads  to  the  variational  principle: 


^|{lpv2  “  0  +  ♦  8  (•^+▼•7#]  +  T'  (~+(v'7)X||<l3rdC  -  0  (1-14) 


where  we  have  included  also  the  additional  constraint  on  the  identity 
of  fluid  particles: 


—  +  (v7)X  »  o 
d  t 


(1-15) 
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or  the  change  in  initial  location  x  of  a  particle  which  is  at  "r"  at 
time  "t"  is  zero  if  one  follows  the  particle’s  motion.  This  may  appear 
tautological  but  is  necessary  to  produce  physical  results  if  one  per¬ 
forms  a  variation  where  P,v  are  fixed  at  the  endpoints  rather  than  the 
13 

displacements.  Without  (1-15)  one  obtains  the  restriction  that 

isentropic  flow  must  be  irrotational. 

Varying  the  independent  quantities  leads  to 

♦  ♦  ♦  *  ♦ 

5v:  pv  -  pV$  +  pVs  +  *  0 


6p: 

6s: 


M  +  v..au_ 

3p  i  av.p 

j 


+  V  -v0  -  0 

at 


3! 

at 


v  ■  va  •  o 


(1-16) 


5X ;  ♦  V  •  ( v*  . )  *  0 

at  i 

whereas  varying  <f>,yf8  leads  to  (1-12),  (1-15),  and  (1-13).  Combining 
equations  (1-16)  leads  to  Newton's  law  of  motion  for  the  dispersive 
continuum 


-VP  + 


a  u _ 

a(ap/ar  j_) 


where  we  have  set 


(1-17) 


au 


P  •  -  0  +  PU  +  psT  +  V.p1 

1  av.p 


,  (1-18) 
and  Introduced  as  an  extension  of  the  1st  or  2nd  laws  of  thermodynamics : 


dU  ■  |A*dp  +  Td(ps)  ♦  -  d(V.p) 

av.p  i 


(1-19) 
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From  these  basic  equations  (1-12,  1-13,  1-17)  we  can  derive  the  laws  of 
momentum  and  energy  conservation 


+  !£ij  -  0  and  ^-<^pv2+U)  +  7-Q  -  0  ,  (1-20) 

at  ar^  at  2 


where 

„  a  au 

p.  .  ■  P6.  .  +  pv.v.  -  p  —  • -  and 

jo  jo  *>3  a^  a(ap/ari> 


Qi 


(p  v  /2  sT)pv .  -  pv. - - 

1  Jar.a?.P 


av.p 

i 


\  »  / 


Note  the  stress  tensor  is  no  longer  diagonal  so  that  Pascal's  law  is 
not  obeyed  by  a  dispersive  system.  For  isentropic  flows  (s»constant) 
Kelvin's  circulation  theorem  and  Lagranges  theorem  follows  as  for  non- 
dispersive  fluids.  For  s  constant  (1-17)  becomes: 

. .t  -  .j-jl) 


so  that 


±-r2  +  p  +  sT  -  7 .  •— — 
2  i  87. p 

l 


is  constant  along  a  streamline  and  for  irrotational  flow: 


7* 


and 


2 

~  +  —  ♦  n  +  sT  -  7.  •— —  ■  constant  (1-21) 
8t  2  1  37.p 


throughout  the  fluid. 


where  Pg  and  ag  are  the  constant  density  and  entropy  of  the  equilibrium 
state,  e  is  a  measure  of  the  deviation  from  equilibrium,  and  the  set 
of  functions  (P^v^  s^)  are  n  order  corrections  to  be  calculated  in  an 
iterative  scheme.  The  fluid  variable  expansions  are  substituted  into 
the  fluid  equations  and  since  the  quantity  £  is  by  hypothesis  small. 
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ve  demand  that  the  equations  be  satisfied  for  each  power  of  £.  Iso¬ 
lating  terms  of  the  first  power  in  e,  we  obtain  the  linearized  fluid 
equations  from  which  the  first  order  solution  (p^.v^.s^)  is  determined 
subject  to  given  initial  conditions.  The  nC^  order  equations  are  linear 

in  the  functions  p  ,  s  and  v  ,  and  are  comprised  only  of  those  P  ,3  , 
mm  ui  00 

and  v  for  which  m  is  less  than  n.  Thus,  once  the  initial  conditions 

HI 

are  specified  we  can  determine  p,s  and  v  to  any  order  of  approximation 
by  iterating. 

The  specific  entropy  is  constant  to  all  orders  of  approximation. 

To  prove  this  assertion  consider  the  n1"*1  order  equation  for  entropy: 


5 

— -  s 
at  n 


v;  *7 


s 

n- 


1 


+ 

+  v 

n 


s 

o 


0 


5 

If  is  constant  for  m  less  than  n  then  s^-0  implying  that  is 

constant.  Since  the  entropy  at  equilibrium  s  is  constant,  s  for  all 

o  n 

n  Is  constant  by  induction.  Thus  the  dependence  on  "s"  will  be  sup¬ 
pressed  from  now  on.  Equations  valid  up  to  second  order  can  be  obtained 

2 

be  keeping  terms  of  order  (5p)  .  Expanding  the  pressure  (1-26)  we 
obtain 


P 


lap 


28) 


Similarly  for  the  energy  density  we  have 


?a  a  Up 

ar. a V.p 
j  1 


[Sl&P.  1 

♦  Sell 

[dhp  ]6p 

♦  P  21 

if 

i56p 

’Wr/ 

ap  1  e 

5ri.arj 

edpl 

evariJ 

drj 

(1-29) 


so  that  (1-22)  becomes  (valid  to  second  order  in  5o) 
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D<Sp)  -SrsT^U^Aj  ‘’"i 

i  i  °p 


ifSa  ] (76p)3  -  *>fM|  jli£ —  .  (p«  ]«£»£.  }  (1.30) 

2^3p  e'  l8p  le'dr^ar;  dp  e'dr^ar;' 


where: 


0  5  2-2  +  cV  ♦  a74  ,  a  =  (p?)  ,  and  c2  5  ]  . 

at  *  « 

The  underlined  terms  in  (1-30)  are  nonlinear  as  well  as  dispersive. 
These  will  be  small  for  the  cases  which  we  consider  as  we  will  be 
concerned  with  the  leading  effects  of  nonlinearity  and  dispersion. 
Thus  the  equation  which  forms  the  main  focus  of  this  thesis  is 


0<6p) 


! - I  ]6.  .  ♦  p  v.v.  ) 

•i&rj  \  2  lap2  |.J  ij  e  i  j  / 


(1-31) 


The  dispersion  arises  from  a  (the  dependence  of  U  on  (^p)2)  and  the 
nonlinearities  are  due  to  the  doppler  shift  and  the  deviation  from 
Hooke's  law.  Employing  the  expansion  (1-27)  the  equation  of  motion  for 


P^  is: 


QPl  *  0 


(1-32) 


from  which  v^  can  be  found  by 


— 1  *  p  7  *v  ■  0 
3t  e  1  u  • 


(1-33) 


The  dispersion  law  is  from  (1-32): 


2  2,2 

<*}  ■  c  k  ♦  ak 


(1-34) 
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into  two  pieces:  one  traveling  at  the  long  wavelength  speed  of  sound 
and  the  other  at  the  group  velocity  of  the  wave  packet.  We  calculate 
the  momentum,  mass,  and  energy  in  these  pulses  and  compare  it  to  the 
properties  of  fluid  solitons  which  are  quite  similar  except  for  sign 
differences.  The  following  table  summarizes  many  of  the  properties  of 
the  dispersive  fluid  pulses  and  solitons. 


In  the  table,  v^  and  u  are  the  group  velocity  and  speed  of  the  packet 

and  soliton,  and  ,  p  de 

i  ■  l  ♦  r — 
e  ie 
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In  Appendix  A  we  also  discuss  Che  relationship  of  Che  pulse  momentum  Co 
Che  Rayleigh  and  Langevln  radiation  pressures. 
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CHAPTER  II 


SCATTERING  OF  SOUND  BY  SOUND  IN  A  DISPERSIVE  CONTINUUM 

la  this  chapter  we  calculate  the  scattering  of  sound  by  a  back¬ 
ground  of  sound  and  show  how  these  effects  are  sensitive  to  the  dis¬ 
persion  of  the  medium.  We  will  use  the  dispersive  fluid  equations  to 
third  order  in  the  amplitude  which  involves  therefore  two  iterations 
of  the  perturbation  scheme.  We  imagine  that  the  background  is  made  up 
of  a  superposition  of  traveling  waves  and  focus  our  attention  on  a 
given  impressed  wave  which  interacts  with  each  of  the  background  waves 
via  the  nonlinearities.  At  the  first  order  of  interaction  the  non- 
linearities  generate  sum  and  difference  frequencies  of  the  Impressed 
wave  and  background.  At  the  next  order  of  interaction  there  is 
generated  waves  at  the  same  wavenumber  and  frequency  as  the  original 
wave.  This  constitutes  a  correction  to  the  original  wave  due  to  its 
nonlinear  interaction  with  the  background  and  yields  the  attenuation. 
For  the  impressed  wave  we  take  for  the  velocity  and  density: 

♦  A  .*►  k  ♦ 

va  *  lcaA ' exP( ika •  r-iu^t )  +  c.c.  and  5pfl  m  .  (2-1) 

where  c.c.  denotes  complex  conjugate  and  uj  (k  )  is  determined  by 

3  3 

(1-34).  Similarly  for  the  background  we  set 

♦  *  *  ♦  k  ♦  ♦ 

v(k)  ■  kB(k) ’exp(ik’r-iut)  c.c.  and  5p(k)  ‘-'P^'k  .  (2-2) 

The  first  order  equations  are  solved  by  the  sum  since  superposition 
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applies  in  Che  linear  approximation: 


v(k) 


♦  v 


and 


P,  ■  5p(k)  +  5p 

1  a 


As  we  are  interested  in  the  change  of  due  to  the  background  we  will 

2 

neglect  terms  which  go  as  (<$p  ) “  which  are  due  to  the  self  Interaction. 

a 

2 

Our  assumption  that  the  dispersion  is  not  too  strong  or  ak  /  -><<1 
causes  the  significant  contribution  to  the  wave  scattering  to  arise 
from  nearly  collinear  waves. 

To  find  the  first  nonlinear  correction  (p^.v^)  we  solve  the 
second  order  equation  which  is  written  below  with  the  first  order 
solutions  (2-1,  2-2)  substituted: 


where 


□  p2  ■  -2Ap  A b£  (k  i  k)  exp( LQ  )  +•  c.c. 
e  ±  a  ± 


4 

+  a  7 


A 


$ 


and  Q  ■  (k  t  k)  -r  -  (u>  ±  u)t 
±  a  a 

The  particular  solutions  for  (Pj,^)  are  found  by  crying  solutions  of 
the  form: 


P2  *  l  V  exp(ifl  )  >  c.c.  and  v2  -  l  W.exp(i Q)  +  c.c. 
+  *  *  XX 


and  then  solving  for  V±  and  W±,  The  homogeneous  solution,  which  is  a 
solution  to  the  linearized  equations,  is  chosen  so  that  the  whole 
solution  satisfies  the  initial  conditions:  the  conditions  being  that 
Pn  for  n>l  vanish  initially  since  our  first  order  solution  is  con¬ 
structed  to  give  the  initial  density  and  velocity  fields.  The  whole 
solution  is  found  to  be: 
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P2  *  Pe2  Ota±  k)&4  +  c.c.  and  v2  -  cl  (ka±  k)R±  +  c.c.  ,  (2-3) 


where 

^  ^  ^  (exp(-i(w  ±w)t)-exp(-iaj±t)) 

R  *  -2AAB(k  ±  k)exp(i(k  ±k)*r)  - - - 5 - = - 

4  a  a  u)+2  -(cu  tu))2 


(ka±  k)  *  ♦ 

w  m  — “ - u(k  ±  k) 

*  lka±  kl  a 


That  the  sound  wave  amplitude  must  change  in  the  process  of  pro¬ 
ducing  (^2*^2^  can  be  Proven  usin8  energy  conservation.  We  find  the 
quantitative  change  from  0^  which  is  determined  by  solving  the  third 
order,  perturbation  equation: 


3  8r 


*  (’.(''UV2j  *  *2iV}  (2-‘’ 


where  we  have  ignored  terms  cubic  in  first  order  quantities  since  their 
effect,  which  is  studied  in  Appendix  B,  is  to  alter  che  sound  wave  speed 
but  not  the  amplitude.  Also  in  Appendix  B  we  show  that  there  are  dis¬ 
persion  effects  that  have  a  logarithmic  dependence  on  the  sound  fre¬ 
quency.  The  attenuation  effects  are  stronger  since  they  result  from 
resonant  mode  coupling  at  both  orders  of  iteration. 

The  change  in  che  sound  wave  amplitude  is  identified  as  the 
amplitude  of  that  harmonic  component  of  p ^  which  has  the  same  wave- 
number,  frequency,  and  phase  as  the  sound  wave.  This  harmonic  component 
is  generated  only  by  those  terms  in  che  differential  equation  corre¬ 
sponding  to  these  same  characteristics.  Therefore,  when  we  write  the 
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equation  (2-4)  with  first  and  second  order  solutions  (2-1,  2-2,  2-3) 
substituted,  we  remove  all  but  the  relevant  terms: 


2  2  t 


*  *"  *  'r'->  r»  22 

P3  ■  -A|B  |4pgckaA  expdlt^-r-iu^t)  ((w±?  u)  -  c.c.  , 


where 


2  cos((<iU--«aH:w)t)  -  1 

Q.  •  (k  i  k)  - - - = - r— 

( ( "(wa-  w)  ) 


The  relevant  solution,  which  is  found  in  the  same  way  as  for  p^,  is 

2  2  2  +  ♦  - 
P3  ■  AlB  Up^ck^A  exp(  ik^’r-i^t)!  Q+  *  c.c. 

The  change  in  the  sound  wave  amplitude  is  the  amplitude  of  this 
solution: 


2  2  2  2  r  7 

4AIB  Ic  A  \  [  Q_  +  Q+  J 


(2-5) 


At  resonance  where  u)±-(uj  ±  a))»0,  the  Q±  grows  quadratically  in  time 

a 

whereas  off  resonance  the  Q±  is  bounded  and  oscillatory.  When  disper- 
ak^ 

sion  is  weak  ( — j  <<^-)  •  resonance  occurs  when  the  angle  between  the 
c 

directions  of  propagation  of  the  sound  and  background  waves  satisfies 


2  3a  2  2  3a  2 

9  s  2  (k  k)  or  (9-it)  s  2  (k  -  k) 
c  *  c  3 


from  which  we  infer  that  for  resonance  the  waves  must  be  nearly  col- 
linear  and  a  must  be  nonnegative. 

To  find  the  change  in  the  sound  wave  amplitude  due  to  nonlinear 
interaction  with  a  background  distribution,  we  integrate  the  change  due 
to  the  single  background  wave  treated  as  a  differential  piece.  To 
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characterize  the  background,  we  let  E(k,n)  be  the  background  energy 

density  per  unit  wavenumber  at  k  and  per  unit  solid  angle  in  the  n 

direction.  We  choose  the  amplitude  B  of  the  single  background  wave  so 

that  its  wave  energy  equals  the  background  energy  in  the  interval  dk  at 

k  and  in  the  solid  angle  <£}  at  k.  Since  the  wave  energy  density  to 
1  2 

second  order  is  y  p  |b|  ,  we  choose  B  to  be  given  by 

IB2  I  -  -  E(k,k)-dk-dQ 
Pe 

Using  (2-5),  we  integrate  the  effect  of  this  differential  piece  to  find 
the  change  in  the  sound  wave  amplitude  due  to  the  background  distribu¬ 
tion: 


where 


fdQ  Ea.lO-Q^ 


The  angular  integration  can  be  approximated  for  long  times  because 
the  quantity  Q±  becomes  peaked  about  resonance.  Since  the  main  contri¬ 
bution  to  the  integral  comes  from  a  solid  angle  centered  at  resonance 
and  whose  width  is  inversely  proportional  to  time,  we  can  approximate 
the  integral  by  expanding  the  integrand  in  the  angle  about  resonance: 


h±  -  -L  d6,2iie' 


E(k,k) 

4kac^(ka±  k) 


sin2 (  (-P±92/2  »  Bj.)t/2  ) 
{~?*92/2  -  IU) 


where 

ckk.  in 

P+  ■  ,  -  and  R.  *  — ’klka±  kl  . 

£  lkai  kl  t  2c  a 
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1  2 

Using  che  substitution  y  —  y  P±9  +R±  and  for  large  times  letting  the 
upper  limit  of  integration  extend  to  infinity,  we  are  left  with 


where 


*Q±  "  -t‘a_ 
8P+ 


I 

-R+t/2 


sin2(y) 

2 

y 


(2-6) 


For  times  such  that  | tR+|  >>  t,  the  value  of  I  depends  on  the  sign  of  a: 


if  a  <  0  ,  then  120, 
if  a  >  0  ,  then  122. 


When  |tR+|  <<  ir  or  if  a-0,  then  I«l. 

That  Westervelt's  prediction  for  the  attenuation  of  sound  in 

liquid  helium  is  for  some  cases  roughly  half  the  experimental  value  is 

simply  explained  by  the  sensitivity  of  the  attenuation  upon  the  nature 
14 

of  the  dispersion.  Using  the  result  (2-6),  the  change  in  the  sound 

wave  amplitude  due  to  the  background  is  AA  -  -tAB,  where  B,  which  is 

defined  as  the  attenuation  coefficient,  is  given  by 

2 


@(ka> 


2»,c  \  I. 


(ka+  k)  A 

— - -  E(k 


A  f  U 

,ka)dk  - 


(k«-  k)  a 

- -  E(k,ka)dk 

k  k 


‘a  (ka-  k) 


E(k 


,ka)dk  |  (2-7) 


Treating  thermal  motion  as  the  background,  Westervelt  used  nondisperaive 

hydrodynamics  for  his  prediction  for  which  a  is  zero  and  the  value  of  I 
4 

is  unity.  However  if  a>0,  then  the  value  of  I  can  double. 

The  attenuation  coefficient  (2-7)  has  three  terms  each  of  which 
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corresponds  to  a  specific  process.  The  first  term  is  the  attenuation 
due  to  those  nonlinear  interactions  for  which  the  sound  wave  and  the 
background  waves  interact  to  produce  waves  at  the  sum  wavenumber  (k+k&) . 
The  second  and  third  terms  correspond  to  the  processes  that  produce 
waves  at  the  difference  wavenumber  (k  -k) .  The  second  term,  which  is 
the  contribution  from  those  background  waves  for  which  ka>k,  is  negative 
and  therefore  the  sound  wave  is  amplified  in  these  processes. 

Media  for  which  a>0  have  a  dispersive  nature  that  has  been  labeled 

2 

anomalous.  Our  experience  with  most  materials  is  that  — —  <  0,  however 

dk2 

upon  examination  of  the  dispersion  relation  (1-34)  it  is  evident  that  a 

,2 

positive  a  implies  that  — y  >  0. 

dk  5 

The  specific  heat  measurements  of  Greywall  indicate  that  the  dis¬ 
persion  of  liquid  helium  at  low  pressures  is  anomalous  for  long  wave¬ 
lengths  . 

With  the  implications  of  the  specific  heat  measurements  at  high 
pressures,^  we  can  explain  by  using  our  calculations  why  the  attenuation 
of  sound  in  liquid  helium  decreases  dramatically  when  the  pressure  is 
raised  past  18  atmospheres.  The  specific  heat  measurements  show  that 
the  dispersion  of  liquid  helium  becomes  normal  (a<0)  at  pressures  above 
18  atmospheres.  However  if  a  changes  from  positive  to  negative,  our 
calculations  have  shown  that  the  value  of  I,  upon  which  the  predicted 
attenuation  is  proportional,  drops  approximately  from  two  to  zero. 

For  long  times  the  decay  of  the  sound  wave  in  the  presence  of  a 
background  is  exponential.  Though  in  our  perturbation  scheme  we 
implicitly  assumed  that  the  sound  amplitude  was  nearly  constant,  we 
extend  its  validity  by  taking  the  amplitude  to  be  time  dependent.  Then 
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in  Che  interval  (t,t+dt)  the  change  in  the  amplitude  is  given  by  the 
differential  equation: 

dA  -  -0-A(t)dt  , 


•St 

which  is  integrated  to  give  A(t)  »  A(o)  e 

The  attenuation  of  sound  in  liquid  helium  at  low  temperatures  is 
understood  as  resulting  from  the  nonlinear  interaction  of  the  sound 
with  the  background  of  thermal  motion.  It  is  known  from  statistical 
mechanics  that  Che  thermal  background  is  isotropic  with  the  energy 
density  per  unit  wavenumber  given  by  the  Planck  distribution. 

To  obtain  a  numerical  value  for  the  attenuation  of  sound  in  liquid 
helium  at  low  temperatures,  we  will  insert  into  our  calculated  attenua¬ 
tion  the  values  for  the  speed  of  sound  (c) ,  the  density  (p)  and  the 
Gruneisan  constant  (A)  which  are  given  by  the  experiments  of  Abraham 
et  al.^  The  predicted  attenuation  will  then  be  compared  to  the  experi 
mental  results**  for  which  the  attenuation  was  measured  for  sound  fre¬ 
quencies  between  10  and  260  megahertz,  temperatures  between  .1®  and  1®K 
and  pressures  between  0  and  25  atmospheres. 

For  the  range  of  frequencies  and  temperatures  covered  by  the  ex¬ 
periments,  the  expression  for  the  calculated  attenuation  Cakes  the  fol¬ 
lowing  approximate  form: 
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(2-8) 


where  k^  is  Boltzmann's  constant  and  h  is  Planck's  constant.  The 

simplification  arises  because  in  the  experiments  the  sound  frequency 

hu) 

and  the  temperature  satisfy  the  condition:  k^T  <<  1.  This  condition 
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implies  that  most  of  the  energy  in  the  Planck  distribution  resides  in 
wavenumbers  much  greater  than  the  sound  wavenumber,  which  allows  us  to 
take  as  zero  the  value  of  the  sound  wavenumber  as  it  appears  in  the 
limits  of  integration  in  the  attenuation  expression  (2-7) .  The  attenua¬ 
tion  then  becomes  proportional  to  the  total  energy  which  is  given  by 
statistical  mechanics. 

Using  the  experimentally  determined  values  of  the  bracketed 
quantities  in  (2-8),  we  find  the  predicted  attenuation  coefficient  for 
sound  in  liquid  helium  at  zero  pressure  as  a  function  of  frequency  (in 
radians  per  second)  and  temperature  (in  degrees  Kelvin)  to  be 
-6  4 

8  ■  1.5*10  Iw  T  (dB/cm) 
r  & 

Figure  1  is  the  plot  of  Che  predicted  (1-1  and  1*2)  and  measured 
attenuations  versus  temperature  for  a  15  megahertz  sound  wave  in  liquid 
helium  at  zero  atmospheres.  (The  figures  are  on  pages  29-32.) 

Figure  2  is  a  plot  of  the  measured  attenuation  versus  frequency  at 
zero  atmospheres  and  .2°  Kelvin.  The  solid  lines  are  the  predictions 
of  our  calculation  (1*2)  and  Westervelt's  (1*1). 

Figure  3  displays  the  measured  and  predicted  (1*1  and  1*2)  attenu¬ 
ations  versus  pressure  for  a  256  megahertz  sound  wave  in  helium  at  .22° 
Kelvin. 

Figure  U  shows  the  measured  attenuation  versus  temperature  for  a 
15  megahertz  sound  wave  at  a  pressure  of  25  atmospheres.  The  predicted 
attenuation  is  zero  because  the  dispersion  is  normal  (ot<0)  above  18 
atmospheres.  For  temperatures  above  .6®  Kelvin,  the  attenuation  is  not 
zero  as  predicted  because  the  sound  is  interacting  with  "rotons"  which 
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are  high  frequency  thermal  components  whose  energy  grows  rapidly  with 
temperature. 

We  can  explain  by  considering  bandwidth  effects  why  the  experi¬ 
mental  points  in  Figures  1  and  2  might  lie  between  our  prediction  (I»2) 
and  Westervelt's  (1*1).  Bandwidth  effects  result  from  the  uncertainty 
in  frequency  inherent  in  any  decaying  waves.  In  our  case  both  the 
sound  and  thermal  waves  are  decaying  and  if  these  decays  are  large  then 
the  approximations  made  in  evaluating  the  integrals  over  the  background 
distribution  become  poor.  In  Appendix  B,  we  use  a  Green's  function 
technique  to  show  that  the  effect  of  this  bandwidth  is  to  multiply  our 
calculated  attenuation  by  the  factor: 


where  T  is  the  lifetime  of  a  typical  background  wave.  Including  this 
factor,  we  find  our  attenuation  is  equal  to  Westervelt's  at  frequencies 

fV) 

where  a  u)aT<<^  an<*  is  twice  Westervelt's  at  sound  frequencies  where 
the  inequality  is  reversed. 

The  frequency  dependence  of  the  attenuation  predicted  by  using  the 
bandwidth  effect  does  not  agree  with  the  experimental  points  of  Figure 
2.  Possible  explanations  for  the  discrepancy  between  theory  and  experi¬ 
ment  may  arise  from  considering  the  following:  misinterpretation  of  the 
data,  long  thermal  lifetimes,  or  attenuation  due  to  interaction  with 
second  sound.  For  the  first  point,  the  experimental  data  may  have 
systematic  error  because  to  find  absolute  attenuations,  risky  assumptions 
were  made  about  the  temperature  dependence  of  the  attenuation? 

Secondly,  it  may  be  Incorrect  to  assume  as  we  did  in  calculating 
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the  bandwidth  effect  that  the  thermal  lifetime  is  shorter  than  the  sound 

wave  lifetime.  Although  the  dispersion  in  liquid  helium  is  anomalous  at 

low  wavenumbers,  it  is  always  normal  at  sufficiently  high  k.  Thus  there 

fd2 

is  an  inflection  point  — ~  -0  at  some  finite  k*k. .  If  the  background 

’■dk 

is  concentrated  at  wavenumbers  greater  than  k^,  then  the  background 
would  be  quite  long  lived,  perhaps  with  a  lifetime  longer  chan  Che 
impressed  sound  wave.  This  case  would  correspond  to  a  dispersion  rela- 
tion  such  as  w  »c  4k4+cdc  -a '  k° .  Although  the  results  will  be  different 
from  those  calculated  here  we  believe  that  the  same  procedure  will  work 
and  may  account  for  the  discrepancy. 

Finally  there  can  appear  in  the  distribution  of  sound  waves  com- 
pressional  modes  of  the  local  energy  density  or  so-called  second  sound. 
The  appearance  of  the  new  mode  is  a  spontaneous  symmetry  breaking  in 
the  nonlinear  hydrodynamics.  Now  there  can  appear  new  nonlinear  pro¬ 
cesses  in  which  the  sound  wave  interacts  with  second  sound  and  that 
can  still  further  complicate  the  comparison  of  these  results  with  the 
properties  of  liquid  helium. 16 
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FIGURE  1.  Attenuation  ve.  temperature  for  o  15  MHz  oound  wave  at 
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PRESSURE  (ATM) 

FIGURE  3.  Attenuation  vs.  pressure  for  a  256  MHz  sound  wave 
at  a  temperature  of  0. 22  K. 


CHAPTER  III 


BOLTZMANN  EQUATION  FOR  SOUND 

In  Che  last  chapter,  where  we  considered  the  nonlinear  interac¬ 
tion  between  a  sound  wave  and  a  background,  we  assumed  that  the  back¬ 
ground  energy  distribution  was  constant  in  time.  In  this  chapter,  we 
determine  the  nonlinear  self- interaction  of  an  isotropic  background, 
and  express  how  the  background  distribution  changes  in  time  via  a 
kinetic  equation  derived  from  the  hydrodynamics.  If  we  impose  a  Debye 
type  cutoff,  the  stationary  background  distribution  of  the  kinetic 
equation  is  a  "white  noise"  distribution.  Next  we  present  a  reasonable 
method  of  adding  the  zero  point  motion  contributions  to  the  kinetic 
equation.  The  form  of  the  zero  point  distribution  is  determined  using 
the  kinetic  equation  and  a  law  analogous  to  Wien's  law  for  radiation. 
With  this  form,  the  improved  kinetic  equation  has  a  stationary  solu¬ 
tion  that  is  a  Planck  distribution.  The  chapter  ends  with  the  presen¬ 
tation  of  a  classical  H-theorem  and  a  discussion  of  some  effects  of 
zero  point  motion. 

To  simplify  the  discussion  we  will  impose  the  restrictions  that 
the  medium  be  nondispersive  and  that  the  background  be  low  enough  in 
energy  so  chat  its  3elf-interaction  can  be  considered  as  being  the 
collection  of  binary  interactions  between  the  harmonic  components  that 
constitute  the  background. 

Our  immediate  goal  is  to  derive  a  kinetic  equation  which  gives 
the  time  rate  of  change  of  u(k),  the  background  energy  per  unit 


wavenumber  ac  k  or  so-called  spectral  intensity.  This  background 
energy  at  k  changes  not  only  by  its  direct  interaction  with  the  rest 
of  the  background  but  also  by  restituting  processes  resulting  from  the 
binary  interaction  of  other  background  pieces.  The  change  due  to 
direct  interaction  is  determined  directly  from  the  attenuation  of  sound 
calculated  in  the  last  chapter.  The  effect  of  the  restituting  pro¬ 
cesses  is  found  with  the  additional  use  of  energy  conservation  which  is 
a  property  of  the  nonlinear  continuum  mechanics. 

For  an  isocropic  background,  che  time  rate  of  change  of  the  back¬ 
ground  energy  at  k  due  to  direct  interaction  is  the  same  as  that  of  a 
sound  wave  of  wavenumber  k  and  equal  energy.  Having  shown  that  for  a 
sound  wave  the  rate  of  change  of  the  density  amplitude.  A,  is  given  by 

tr  “V  ■  -  s<k.>*‘v  . 

and  knowing  that  energy  is  quadratic  in  amplitude,  we  deduce  that  the 
time  rate  of  change  of  the  background  energy  density  u(k)  for  direct 
interaction  is  expressed  by 

—  «(k)  -  -2{3(k)u(k) 
dc 

where  B(k)  in  terms  of  u(k)  was  derived  in  the  last  chapter: 

"  8p~cA  (I  ♦  |  u( q)dq  -  j  u(q)dqj. 

To  complete  the  kinetic  equation,  we  must  find  the  rate  of  change 
of  the  background  energy  density  due  to  restitution.  When  two  waves 
interact  they  produce  waves  whose  wavenumbers  are  equal  to  the  sum  and 
difference  of  the  wavenumbers  of  the  interacting  waves.  It  follows 
that  in  a  restituting  process,  the  background  energy  at  k  changes 
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because  Che  wavenumber  k  is  either  the  sum  or  difference  of  Che  wave¬ 


numbers  of  many  pairs  of  directly  interacting  background  pieces.  The 
changes  due  to  the  sum  and  difference  cases  will  be  treated  independ¬ 
ently  and  will  be  calculated  for  single  pairs  of  interacting  background 
pieces  by  using  energy  conservation  from  which  we  know  that  the  energy 
gained  in  restitution  is  the  energy  lost  by  direct  interaction.  The 
total  restitution  will  be  given  by  an  integral  over  all  possible  in¬ 
teracting  pairs. 

We  begin  by  focusing  on  the  single  sum  process  in  which  the  back¬ 
ground  energy  in  the  interval  AK  at  wavenumber  (k-q)  Interacts  with  the 
energy  in  the  interval  dq  at  q  to  restitute  energy  to  the  background  at 
the  sum  wavenumber  k.  The  time  rate  of  energy  loss  of  these  interac¬ 
ting  waves  can  be  obtained  by  dissecting  our  kinetic  equation  for 
direct  interaction,  which  is 

~  «<k)  *  -20(k)u(k)  .  (3-1) 

a  t 

The  attenuation  coefficient  for  a  sound  wave  of  wavenumber  k  resulting 
from  just  sum  interaction  with  a  background  is  given  by  the  first  term 

in  the  brackets  of  the  expression  for  B(K): 

•  2 

0  (K)  *  ^  f  ■  u(q)dq  ,  where  Q  »  — —  A 

sum  2  J0  <1  4pec 

Singling  out  the  contribution  of  the  energy  in  the  interval  dq,  sub¬ 
stituting  into  Eq.  (3-1)  and  letting  K-k-q,  we  find  that  the  time  rate 
of  change  of  the  background  energy  in  A K  at  (k-q)  for  our  single  sum 
process  is 

—  u(k-q)AK  *  -<}  —  u(q)dq  u(k-q)AK 
dt  q 
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Likewise,  the  rate  of  change  of  the  energy  In  dq  at  q  is 

2 

—  u(q)dq  ■  -Q  k  ■  11  u(k-q)AR  u(q)dq 

dt  (k  -  q) 

Since  the  energy  lost  by  both  waves  In  direct  sum  Interaction  is  gained 

In  restitution  in  the  interval  AK  at  k,  we  deduce  that 

2 

—  u(k)  ■  Q  -  k - -  u(k-q)u(q)dq  . 

dt  q(k  -  q) 

Integrating  this  expression  over  all  possible  pairs  of  interacting 
waves,  we  have  the  rate  of  change  of  the  background  energy  at  k  caused 
by  the  restituting  sum  interactions: 

k/2  3 

u(k)  ■  Q  [  - - -  u(k-q)u(q)dq  .  (3-2) 

dt  Jg  q(k  -  q) 

The  contribution  to  the  kinetic  equation  from  the  restituting  dif¬ 
ference  processes,  which  is  calculated  in  the  same  manner  as  the  sum 
processes,  is 

•  3 

•i-  u(k)  »  Q  f  - - -  u(k*q)u(q)dq  .  (3-3) 

dt  ^0  qvk  ♦  q) 

In  calculating  this  contribution,  we  use  the  attenuation  coefficient  for 
sound  interacting  with  a  background  via  difference  processes: 

0dif(k)  *  j  {  I0  'k  q  ^  u(q)dq  "  fk  •  (3-4) 

The  significance  of  the  two  terms  is  that  when  two  waves  interact  in  a 
difference  process,  the  wave  of  lesser  wavenumber  gains  energy  while 
the  wave  of  greater  wavenumber  loses. 

The  complete  kinetic  equation  obtained  by  adding  the  contributions 
of  the  direct  process  (3-1)  to  that  of  the  sum  and  difference 
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restituting  processes  (3-2,  3-3)  is 


—  —  u( k)  »  /  -£(k)u(k)  ♦  [ 

Q  dt  V  - - -  1, 


k/2  3 

k 


o  q(>«  -  q) 

3 


u(k-q)u(q)dq  + 


f  __Ji - _  u( k+q)u(q)dq  )  (3-5) 

Jo  q<k  ♦  q)  > 


We  will  refrain  from  application  of  this  equation  until  ve  improve  it 

with  the  inclusion  of  zero  point  motion  contributions.  However,  we 

mention  chat  if  a  wavenumber  cutoff  is  imposed  analogous  to  that  which 

Debye  argued  for  solids,  then  Che  kinetic  equation  has  a  stationary 

distribution  (i.e.,  -0)  of  the  form: 

dt 

2 

u(k)  ■  constant  *  k  , 


which  is  representative  of  the  statistical  mechanic  notion  of  equiparti- 
tion  of  energy. 

Our  attention  is  now  directed  toward  determining  how  the  kinetic 
equation  is  altered  by  zero  point  motion.  Though  the  effects  of  zero 
point  motion  has  traditionally  been  thought  to  be  outside  the  domain 
of  classical  physics,  we  will  show  that  it  is  possible  to  introduce 
2ero  point  motion  in  a  classical  vein.  The  zero  point  motion  is  postu¬ 
lated  as  an  immutable  energy  distributed  uniquely  over  all  wavenumbers, 
existing  even  in  the  absence  of  thermal  motion,  but  nevertheless  able 
to  interact  in  some  restricted  sense  with  other  waves. 

The  implementation  of  our  postulate  becomes  apparent  when  we  con¬ 
sider  a  sound  wave  in  the  presence  of  a  zero  point  wave.  We  would  like 
to  treat  the  interaction  of  the  sound  wave  with  the  zero  point  wave  in 
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the  same  way  as  any  two  waves,  but  since  the  energy  In  a  channel  which 
has  only  zero  point  energy  cannot  diminish,  we  must  somehow  restrict 
the  interaction.  For  the  case  of  two  waves  interacting  in  a  sum  pro¬ 
cess,  both  waves  lose  energy;  therefore  if  one  of  the  waves  is  a  zero 
point  wave  whose  energy  cannot  diminish,  the  sum  process  must  be  elimi¬ 
nated.  In  the  difference  process,  the  wave  of  lesser  wavenumber  gains 
energy  and  the  wave  of  greater  wavenumber  loses.  As  in  the  sum  process, 
the  possibility  that  a  wave  of  greater  wavenumber  in  the  difference 
process  be  a  zero  point  wave  must  be  eliminated,  but  we  may  allow  a 
sound  wave  to  interact  with  a  zero  point  wave  of  lesser  wavenumber  if 
Che  energy  normally  gained  at  the  zero  point  wavenumber  is  placed  into 
sound  wave  energy  at  that  wavenumber. 

To  calculate  the  zero  point  contributions  to  the  kinetic  equation, 
we  utilize  that  the  interactions  of  the  background  waves  with  the  zero 
point  waves  is  just  a  subset  of  the  interactions  of  the  background 
self-interaction  whose  contributions  to  the  kinetic  equation  have  been 
calculated.  The  time  rate  of  change  of  the  background  energy  at  k 
caused  by  zero  point  motion  is  due  to  difference  processes  of  three 
types: 

1.  The  direct  Interaction  of  the  background  energy  at  k  with 
zero  point  waves  of  wavenumber  less  chan  k, 

2.  The  direct  interaction  of  the  zero  point  energy  at  k  with 
background  waves  of  wavenumber  greater  than  k, 

3.  Restitution  from  the  background  energy  at  (k+q)  interacting 
with  the  zero  point  energy  at  q  (for  all  q). 

For  each  of  these  types,  the  time  rate  of  change  of  the  background 
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energy  ac  k  caused  by  the  background  has  been  calculated  to  be,  respec¬ 
tively, 

k  2 

(1)  -  u(k)  f  u(q)dq  , 

iQ  « 

•  2 

(2)  u(k)  f  1! u(q)dq  , 

jk  q 

•  3 

(3)  f  - - — — ■  u(k+q)u(q)dq 

J0  q(k  +  q) 

As  an  example  of  how  this  procedure  works  consider  the  change  in 
energy  in  channel  k  due  to  sound  in  channel  q<k  due  to  difference  pro¬ 
cesses: 

~  u(k)dk  -  -Q  u(  q)dq  u(k)dk  ,  (3-6) 

where  we  have  used  (3-1,  3-4) .  Similarly  for  q>k,  we  have 

2 

•i-  u(q)dq  *  Q  -- -  u(k)dk  u(q)dq  .  (3-7) 

dt  k 

Next  on  the  right  hand  sides  of  th«.se  expressions  we  replace  u  with  u+v 
and  focus  attention  on  the  contribution  linear  in  the  zero  point  energy 
per  unit  wavenumber,  which  we  will  label  v,  (the  terms  quadratic  in  u 
have  already  been  incorporated  into  (3-5)): 

2 

~  u(k>  •  -Q  -  [  u(q)v(k)  ♦  u(k)v(q)  ]  dq  (3-8a) 

2 

u(q )  *  Q  [  u(q)v(k)  +  u(k)v(q)  ]  dk  (3-8b) 

dt  k 
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so  that  from  (3-8a)  we  see  that  the  process  whereby  a  zero  point  fluc¬ 
tuation  of  wavenumber  k  interacts  with  an  excitation  of  wavenumber  q  is 
forbidden  since  it  removes  energy  from  channel  k.  Deleting  this  term 
from  (3-8a)  and  (3-8b)  yields  direct  and  restituting  contributions  to 
the  change  in  energy  density  given  by  expressions  (3-6)  and  (3-7)  above 
with  u(q)  replaced  with  v(q).  Replacing  also  the  background  energy 
density  of  lesser  wavenumber  in  the  8(k)  in  (3-1)  with  the  zero  point 
energy  density  (i.e.  u(k)-*v(k))  we  find  the  zero  point  contribution  to 
the  kinetic  equation  which  combined  with  (3-5)  yields  the  complete 
kinetic  equation^  ^  d 


u(k)  -  {-u(k)[g(k)  *  [Q  v(q)dq) 


I, 


k/2  3 

k 


q(k-q) 


-  i 

u(k-q)u(q)dq  ^ ■  u(q)u(k+q)dq 

J0  q(k+q) 

3 


♦  v 


(k)[  u(q)dq  +  [  ------  v( q)u(k+q)dq  ). 

Jk  q  ^0  q(k+q)  / 


The  stationary  solution  of  the  kinetic  equation  is  the  background 
equilibrium  distribution.  To  find  this  stationary  solution,  we  must 
first  determine  the  form  of  the  zero  point  distribution. 

Using  the  kinetic  equation  and  a  fluid  law  analogous  to  Wien's  law 
for  radiation,  we  are  able  to  show  using  the  second  law  of  thermody¬ 
namics  that  the  zero  point  distribution  is  of  the  form: 


3 

v(k)  ■  Ack 


where  A  is  an  undetermined  constant.  The  fluid  law,  which  is  proven 


in  Appendix  C,  states  that  the  background  equilibrium  distribution 

must  have  the  form:  3 

u(k)  •  ck  f(ck/T)  , 
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where  the  function  f 


is  as  yet  unknown.  Substituting  this  form  into 
the  kinetic  equation  we  find  that  ail  the  terms  involving  quadratic 
powers  of  u  yield  expressions  which  are  homogeneous  functions  of  T  and 
k  of  the  eighth  degree.  That  is  replacing  T  and  k  with  aT  and  ak  yields 

Q 

a  times  the  original  expression.  Thus  if  there  is  to  be  a  non- trivial 
balance  between  the  zero  poinc  terms  and  the  thermal  terms  the  expres¬ 
sions  involving  v  must  also  be  homogeneous  of  degree  8  which  means  that 
v  is  proportional  to  k^. 

The  equilibrium  background  distribution,  which  is  the  stationary 
solution  to  the  kinetic  equation  subject  to  Wien's  law,  is  a  Planck 
distribution: 

3  _t 

u(k)  -  2Ack  (exp(yck/T)  -  1)  , 

where  y  is  a  second  undetermined  constant.  That  this  is  a  stationary 
solution  to  the  kinetic  equation  can  be  shown  directly  with  the  follow¬ 
ing  identities: 

f(q)f(k+q)  *  f(k) f(q)  t  f(k)f(k+q)  t  f(k+(q±q)/2)  , 

where  f(k)  -  (exp(ak)-l)  1.  Of  the  two  constants,  A  and  Y,  only  one 
must  be  found  by  experiment  since  at  high  temperatures  the  equipartition 
of  energy  determines  the  ratio: 


A  _  kb 
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Einstein  and  Stern  by  considering  the  scattering  of  a  charged 

particle  by  an  electromagnetic  field  attempted  to  derive  the  Planck 

distribution  from  classical  mechanics  by  assuming  a  particular  zero 

point  noise  spectrum.  Their  results  were  off  by  factors  of  two.  Here 

by  considering  self-consistently  the  scattering  of  sound  by  sound  we 

have  arrived  at  the  proper  thermal  distribution  as  well  as  the  zero 

point  spectral  intensity.  These  arguments  have  now  been  extended  to 

19 

the  nonlinear  electromagnetic  equations  by  Roberts. 

That  an  "H-theorem"  exists  can  be  proven  with  the  kinetic  equation. 
If  we  let  H  be  defined  as 


H  **  f  d( lnk)(  (u(k)+2v(k) ) * ln( u(k)+2v( k) )  -  u(k) • ln(u(k) )  ]  , 

0 


then  using  the  kinetic  equation  and  the  cubic  form  of  the  zero  point 
motion  we  can  show  that 


dH 

dt 


iq'M  "«  [  *d,  “(<!>»(»»<»-))  u-b, !„(,/„) 

J0  A  qK(K-q) 


3 

and  b  -  (1  +  K  /u(K)) 

Since  the  integrand  is  always  positive,  we  conclude  chat  4^-  >  0,  which 

at 

is  the  H  theorem.  Demanding  that  H  is  stationary  and  imposing  Wien's 
law  for  fluids,  we  again  obtain  the  background  equilibrium  distribution. 

As  a  final  application  of  the  kinetic  equation,  we  calculate  the 
attenuation  and  decay  spectrum  of  a  sound  wave  interacting  with  zero 


where 


•  ( 


1  + 


u(q) 


3  .]-(i  ♦  iEaL.j 

‘  K  u(K-q)j 


point  motion, 
follows: 


The  sound  wave  is  represented  in  terms  of  u(k)  as 
u(q)  ■  ug6(kg-  q)  , 
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where  uq  is  the  sound  wave  energy  and  kQ  is  the  sound  wavenumber.  Sub¬ 
stituting  this  u(q)  into  the  right  hand  side  of  the  kinetic  equation, 
we  find  that  the  first  zero  point  term  gives  the  attenuation  of  the 
sound  wave: 


dua 

dt 


-  «o 


30 


5 

*9 


A2h  5 
480p^ir^  ® 


and  the  last  two  zero  point  terms  give  the  spectrum  of  decay: 


duo  /  2  3 

T3*  "  uqQAc  (kg  -  k)  k  / kg  for  k<kn  and  zero  for  k>k„. 
at  U  U 


Landau  and  Lifshitz  have  suggested  that  zero  point  motion  effects 
for  macroscopic  variables  can  be  calculated  by  simply  adding  the  zero 
point  energy  to  thermal  energy.  However,  this  method  would  yield  an 
infinite  sound  attenuation  in  contrast  with  the  method  presented  in 
this  chapter  for  which  the  results  agree  with  quantum  theoretical 
calculations. 

The  spontaneous  decay  of  sound  was  measured  for  the  first  time  in 
20 

1981  in  an  experiment  in  which  high  frequency  sound  was  generated  in 

doped  calcium  fluoride  through  the  use  of  laser  pulses  to  stimulate 

non-radiative  transitions.  The  sound  waves  were  detected  through  the 

fluorescence  they  induced  in  the  stress-split  energy  levels  of  the 
2+ 

Eu  doping  ions.  Phonon  lifetimes  were  obtained  from  the  time  de¬ 
pendence  of  the  fluorescence  radiation  and  followed  the  k5  law. 
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APPENDIX  A 


RAYLEIGH,  LANGEVIN  AND  SOLITARY  WAVE 
RADIATION  PRESSURE  IN  A  DISPERSIVE  CONTINUUM 

In  this  appendix  we  use  dispersive  nonlinear  hydrodynamics  to 
determine  the  motion  and  integrals  of  the  motion  for  localized  travel¬ 
ing  disturbances  in  one  dimensional  continuous  media. 

Consider  a  disturbance  which  at  first  order  can  be  expanded  in  a 
Fourier  integral 

m 

Pl(x,t)  -  |  dk-Pl  OOexpdkx-iut)  (A-l) 

with 

w(-k)  ■  -<»>(k)  and  PL(-k)  -  p*(k) 

so  that  the  disturbance  propagates  in  the  "+x"  direction.  Similarly 

m 

J^x.t)  ■  J  dk1  (k)exp(  ikx-iwt) 

— m 

■  f  dk  p, (k)exp( ikx-iwt) 

_•  k  1 

The  inverse  formula  is: 

.  • 

~  p j.  (*> C)exp(-ikx) 
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Substituting  (A-l)  into  (1-36)  yields  for  the  second  order  contribution 
to  the  density 


^2. 

.2  ->* 
2^  P2  +  ^  P; 

dt2 

5x2  Sx4 

II  dkdk'A(k,k'  )exp(  i(k+k'  )x-i(w-K*i'  )t)  ,  (A-2) 


where 

2  <-^A 

A(k,4c')  -  -(kfk'r  P1(k)p1(k') 


and 


A 


1 


£ 

c  dp 


The  solution  for  P,  which  vanishes  at  t»0  and  corresponds  to  a  wave 
moving  in  the  same  direction  as  the  original  pulse  is: 


p„  -  ffdkdk' - A( k , k ' ) exp ( i ( k*k^) xj - _{  exp(-i(umd'  )t)-exp(-iw(k*k'  )t)  } 

*•*  )  +c  (k+k' ;*a(k+k') 


which  is  made  up  of  a  particular  plus  homogeneous  solution  to  (A-2) 


p2  '  P2,p  +  P2  h 


Similarly 


J2  "  J2,p  +  J2,h 


-  [ [ dkdk'-- “y  >«) — 

**  -(v+v')  +c  (k+V)  +a(k+k') 


(  exp^^u^w'  )t)-  ^--1<-,>exp(-iu)(k-t-k'  )t)) 

v  k+k'  (k+k') 


u( k+k’ ) 


where  in  each  case  the  term  in  exp[-i(urhi>')t]  is  the  particular  solution. 
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The  distinction  between  homogeneous  and  particular  solution  has  a 
physical  as  well  as  matheisatical  importance.  The  frequency  and  wave 
number  of  the  particular  solution  are  not  related  by  ohe  dispersion 
relation  (1-34).  Therefore  the  particular  solution  is  a  property  of 
the  traveling  first  order  wave.  It  is  nonzero  only  in  those  regions 
where  the  first  order  disturbance  is  nonzero.  The  homogeneous  solution 
however  is  an  actual  sound  wave  that  propagates  at  the  phase  velocity 
u)(k+k' ) / (k+k' )  and  can  therefore  separate  (for  dispersive  systems)  from 
the  original  first  order  disturbance. 


By  choice  of  the  initial  condition 

jpjdx  ■  o 

but 


(A-4) 


jj2dx  -  2»A  ^-Jdk  |p2(k)|  (A-5) 

fp2  pdx  -  2tA  fjdk  |p2(k) j+(v2(k>  -  c2>  (A-6) 

fJ2,pd*  "  2*AfJdk  lpJ(k)Ku)*(vg(k>  "  (A-7) 


where  the  group  velocity  v^-dw/dk.  Expanding  the  energy  density  (1-23) 
to  second  order: 

[  0  +  jpv2  ]2  -  nep2  *  (c2/2pe)pj  ♦■|pev2  ^i9e(7Pl)2  (A-8) 
leads  to  an  energy  due  to  the  sound  field: 
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Jdx[  U  +  -|pv2  ]2  -  njdx-p2  +  fdx-Pevf 

■  n^Jdx'pj  ♦  -2-p-  Jdk|p2| (l+ak2/c2)  .  (A-9) 

From  (A-4,5,9)  we  find  to  leading  order  (neglecting  corrections  of  order 
ak2/c2)  that  the  ratio  of  total  energy  to  total  momentum  of  the  sound 
field  is 

2c 

o  dc  (A-10) 

1  ♦  -  — 
c  dp 

21 

which  is  the  result  first  presented  by  Rayleigh.  Contrary  to  many 
1  22 

references  *  the  ratio  is  not  c  (except  for  a  monatomic  ideal  gas 

3 

which  has  the  equation  of  state:  PV  -CONSTANT) .  The  following  comment 
should  be  made  about  our  calculation  of  the  energy  and  momentum  of  sound 
waves.  We  have  considered  only  those  effects  which  are  second  order  in 
the  field.  It  is  entirely  possible  that  the  first  order  solutions  p^.J^ 
will  have  a  net  mass,  energy  and  momentum.  However,  we  have  here 
focused  upon  the  quadratic  contributions  since  those  are  intrinsic  pro¬ 
perties  of  the  acoustic  field  when  the  system  is  driven  in  an  oscil¬ 
latory  fashion. 

The  momentum  of  the  acoustic  field  (A-5)  is  proportional  to  A  and 
vanishes  when  the  equation  of  state  is  such  that 

c  •  conatant/p 

For  a  one  dimensional  lattice  of  masses  on  harmonic  springs  we  have 
seen  that  c«/$/n  a  */55n/p  which  provides  the  proper  explanation  of 
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the  statement  that  "phonons  have  no  momentum."  Normally  this  assertion 

is  based  upon  the  reasoning  that  in  a  crystal  lattice  as  a  wave  passes 

each  particle  returns  to  its  initial  location  and  therefore  the  wave 

23 

carries  no  momentum.  However,  the  fact  that  there  is  no  drift  in 
particle  position  merely  requires  /Jdt*0  which  is  quite  different  from 
the  momentum  of  the  wave  which  is  given  by  / Jdx  which  is  proportional 
to  A  and  vanishes,  only  for  the  particular  equation  of  state  appro¬ 
priate  to  the  usually  quantized  harmonic  lattice. 

Next  consider  a  first  order  disturbance  which  is  peaked  so  that 

its  spread  in  wave  number  is  much  less  than  the  wave  number  "k  ".  The 

3 

homogeneous  solution  will  separate  from  the  disturbance  and  its 
particular  solution  so  that  in  time  three  pulses  will  appear :  one  with 
a  mass  and  momentum  given  by  (A-6, 7)  traveling  at  v  ,  another  with  a 
mass  given  by  the  negative  of  (A-6)  and  a  momentum  given  by  the  dif¬ 
ference  of  (A-5)  and  (A-7)  traveling  at  c  and  finally  a  pulse  traveling 

at  v  (2k  )  which  has  zero  mass  and  momentum  to  this  order.  These 
g  o 

assertions  follow  directly  from  substituting  a  peaked  distribution  for 
Px(k). 

Even  when  there  is  no  net  mass  in  the  original  disturbance,  the 

dispersion  separated  pulse  contains  a  net  mass  given  by  (A-6).  Thus 

the  total  energy  of  this  pulse  involves  the  chemical  potential  (A-9) . 

2  2 

Neglecting  corrections  of  order  ak  /c  in  (A-9)  ,  we  find  that  this 
pulse  has  the  total  energy: 


S 
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An  important  distinction  can  be  made  as  regards  the  two  terms  that  con¬ 
tribute  to  the  total  energy.  One  of  the  terms  (the  first  in  the  inte¬ 
gral)  should  be  regarded  as  a  mechanical  energy  since  dissipative  pro¬ 
cesses  can  eventually  convert  this  part  into  heat.  The  other  term  (in¬ 
volving  pg)  will  always  remain  at  this  value  regardless  of  the  irrever¬ 
sible  effects.  Forming  now  the  ratio  of  momentum  to  mechanical  energy 
for  the  separated  pulse  we  obtain 


VgdtQ  )A 

,2  1 

VgOtQ  )  -  c 


( A- 1 2 ) 


We  note  that  when  v  <c,  <0  and  if  v  >c,  0-  > 0  and  the  pulse  moves 

g  2,p  g  2,p 

faster  than  c. 

3 

Bassett  and  Pryce  considered  the  quantum  theory  of  one  dimensional 
dispersive  waves  in  an  anharmonic  lattice.  Though  they  obtained 
Eq.  (A-12)  they  did  not  mention  the  pulse  traveling  at  v^(2kQ)  probably 
because  in  the  quantum  perturbation  scheme  employed  by  them  the  expec¬ 
tation  values  of  energy,  momentum  vanished  for  that  pulse.  Furthermore, 
the  contributions  to  the  total  energy  proportional  to  ve  were  not  in¬ 
cluded  in  their  work.  This  difference  comes  about  because  they  con¬ 
sidered  the  propagation  of  sound  relative  to  an  equilibrium  lattice  that 
is  unstressed.  For  such  a  lattice  Ug-0;  as  can  be  seen  from  the  rela¬ 
tion  (VIZ:  18) 

E  -  UV  «  Mp  -  PV  , 


where  E,V,M,  are  Che  total  energy,  volume,  mass  of  the  system  in  the 
mechanical  equilibrium  state  (^p»0) .  As  is  commonly  done  we  take  E*0 
for  the  mechanical  equilibrium  state.  Thus  if  F  is  also  zero  in  this 
state  we  find  chat  u  must  also  vanish. 
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Had  we  included  the  underlined  terms  in  Eq.  (1-30)  and  so  cal- 

2  2 

culated  corrections  to  the  momentum  (A-5)  of  order  ctk  /c  we  would  have 

o 

found 


(A-13) 


2 

where  of  course  (p^(k))  is  peaked  near  kQ.  The  corrections  proportional 
2 

to  k  still  vanish  for  the  ID  harmonic  lattice  because  in  that  case 
o 

f^a/ ma^l/ p  .  In  general  the  terms  will  not  vanish. 

The  relation  (A-12)  for  the  ratio  of  momentum  to  energy  has 

24 

strong  similarities  to  the  properties  of  solitons  which  are  solutions 
to  Eq.  (1-36)  in  one  dimension  that  have  a  stationary  form: 

p  «  Pe  +  6p(x-ct)  (a-14) 

From  the  continuity  law  this  requires  the  mass  flow  to  be 
J(x-ut)  -  u6p(x-ct) 
so  that  Eq.  (1-36)  yields 


(u2  -  c2)6p 


♦  a 


S26p 

5x2 


(A-13) 


where  terms  cubic  in  <5p  have  been  dropped  along  with  two  space  deriva¬ 
tives.  This  can  be  integrated  directly  to  yield 


where  we  have  taken  a>0  so  that  u<c .  Thus  when  the  sound  pulse  travels 
the  soliton  moves  slower  than  c,  so  that  the  mass  contained  in 
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Che  sound  pulse  has  opposite  sign  to  that  of  the  soliton.  From  (A-16) 
we  see  that  the  small  quantity  for  a  soliton  is  c-u.  If  we  set 

c  -T -  e  <  1  (A-17) 

c 

then 


6p(0,0) 


c  ♦  higher  orders 


in  e 


d26p 

dx2 


i££2e2 

clA 


*  higher  orders  in  e 


(A-18) 


so  that  differentiation  of  3ome  property  of  a  soliton  increases  its 
order  of  smallness.  For  a  sound  wave  both  amplitude  and  phase  are  in¬ 
dependent  but  for  the  soliton  the  two  are  related  by  the  requirement 
(A-14) .  In  view  of  (A-17)  we  were  justified  in  substituting  c  for  u  on 
the  right  hand  side  of  (A-15)  in  order  to  arrive  at  (A-16)  valid  to  this 
order  of  accuracy.  Similarly  to  the  order  of  e  used  here  the  underlined 
terms  in  Eq.  (1-30)  are  higher  order. 

The  momentum  of  the  soliton  is 


uA  £2 

2  2  P 

u*  -  c* 


Jdx(6p)2 


(A-19) 


Introducing  the  mechanical  energy  as  before  we  find  that  the  ratio  of 
momentum  to  mechancial  energy  for  a  soliton  is 


uA 


j2  —  r2 


( A-20 ) 


With  ct>0  a  soliton  exists  with  negative  mass  moving  slower  chan  c. 

There  is  another  wave  packet  which  has  kinematic  properties  very  similar 
to  this  soliton  (see  Eq.  A-12)  which  moves  faster  than  c  and  however 
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spreads  (disperses)  very  slowly  in  time.  The  inverse  scattering 
25 

method  makes  very  general  statements  about  how  initial  disturbances 

evolve  into  solitons.  If  the  initial  pulse  has  the  appropriate  sign 

then  at  a  sufficiently  later  time  (neglecting  dissipation)  at  least  one 

soliton  must  appear.  Thus  one  might  wonder  if  the  dispersion  separated 

pulse  will  eventually  turn  into  a  soliton.  The  answer,  however,  is  no 

for  as  we  have  seen  the  net  mass  in  the  acoustic  pulse  has  opposite 

sign  from  the  soliton  appropriate  to  the  given  sign  of  a. 

Another  quantity  of  interest  in  acoustics  is  the  so-called  Langevin 
26 

radiation  pressure  which  is  the  net  force  exerted  on  a  unit  area  of 
side  wall  as  a  wave  runs  down  the  tube: 

Fl  -  jdt(5P) 

where  to  second  order  the  pressure  is 


6P  -  C2p  +  c  |£.  fi2  <7P,> 

2  ho  *  2  Jn  1 


ap 


dp 


Using  (A-3)  for  0 ^  we  find  for  FL 


'l ■’?.(! 7 •*) 


which  is  the  usual  result.  The  contribution  to  the  Langevin  pressure 
from  the  dispersion  separated  pulse  is  given  to  usual  order  of  accuracy 
by: 


L,P 


2jdt-p2 


+  c 


dc  f  .  2 
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-  It 

[{ ‘‘ 

i  .  ^ 

2cipJ  dp 

+  c  8c  >  p*(k)dk  .  (A-21 ) 

1  Pe 

v  (v2  -  c2) 

g  g 

+  °  dp  1 

2  2 

To  leading  order  in  ak  /c  Che  cerm  in  A  dominates  Che  others.  We  have 

included  Che  term  in  dp^f/d 0  which  comes  from  Che  underlined  term  in 

2  2 

Eq.  (1-30)  since  ic  is  of  Che  same  order  in  ak  /c  as  Che  usual  cerm 
cdc/dp.  For  a  harmonic  lattice  the  cerm  in  cdc/dp  is  the  leading  cerm. 
Ic  muse  be  noted  Chat  although  this  result  becomes  singular  in  the 
limit  v^-*c  it  also  requires  more  and  more  time  for  the  wave  packet  to 
separate  in  this  limit. 

A  connection  can  be  made  between  the  nonlinear  continuum  coef¬ 
ficient  A  and  the  anharmonic  terms  in  the  one  dimension  lactice. 
Generalizing  (1-3)  to: 


r  m  *2  _ 

)  1  *J  *  7 


( A-22) 


where  the  potential  energy  is 


v-I  .1 

n-2  i-l  n> 


i+i  -  *i  -  vn 


(A-23) 


where  x^-position  of  the  particle  which  in  mechanical  equilibrium  is 
located  at  "iao"  where  is  the  unstressed  (F-0)  equilibrium  separa¬ 
tion.  By  definition  the  are  independent  of  "P".  If  a  steady  ex¬ 
ternal  force  "P"  acts  on  this  system  the  equilibrium  spacing  will  change 
to  "a"  such  chat 


P 


-  I 

n-2 


*n 

(n-l)l 


(a  -  aQ) 


n-l 


(A- 24) 
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determines  a(P).  The  speed  of  sound  for  this  one  dimensional  anhar- 
monic  lattice  is 


2  dP 

c*  *  — 

ap 


a2ap  a2  r  *  .  . 

’  ■  3a  »  n.2  ( n-2 ) !  (*  " 


n-2 


(A-25) 


and 


1  |  P  dc  ,  x  _  a  dc 
c  dp  c  da 


“  -  —  I 


*n 


2ac2  n“3  (Q“3)! 


(a  -  a^ ) 


n-3 


(A- 26) 


In  the  stressed  lattice  we  can  refer  the  to  the  separation  "a" 

by  rewriting  (A-23)  as 


V  -  I  \  “  («i+1  -  «i  -  •  ♦  *  -  a0)n 

n-2  i-1  n! 


(A-27 ) 


and  using  the  binomial  expansion  coefficients  c  to  write 

am 


„  7  «B<«>  , 

v "  i  — r  (xi>i 

n-0  n5  1 


xi  -  a  )n 


(A-28) 


where 


r  ®m  ,  .m-n 

~  J.  v 


(A-29) 


yields  the  effective  spring  constants  of  the  stressed  lattice  in  terms 
of  the  spring  constants  of  the  unstressed  lattice.  In  particular 


<L 


®5U)  "  2  l  ~7  cm9(a  "  an> 


m-2 


m-2 


I  m2 


and 


#,(«)  •  6  I  —  c  (a  -  a  )®"3 
3  *-3  ®!  b3  0 
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and  Che  term  in  $  adds  a  constant  to  the  potential  whereas  the  term 
o 

in  <t>i  corresponds  to  a  telescoping  series  which  also  contributes  zero 
to  the  potential.  Note  also  chat: 


frSq  (a) 

da 


4>1  (a) 


3®l(a) 

da 


-  «2(a) 


d®2(a) 

and  -  -  4,  (a) 

da  3 


(A-30) 


Since 


(a)  -  -  P  , 


we  find  for  the  sound  velocity 


c2  -f  •_(.) 

m  i. 


(A-31) 


and  for  Che  Gruneisen  coefficient 


+  £  dc  _  _  a®3<a) 
c  dp 


(A-32) 


2*2(a) 

which  relates  the  hydrodynamic  nonlinear  parameter  to  the  cubic  anhar- 
monic  coefficient  of  the  lattice  for  the  one  dimensional  problem. 

A  precise  transition  from  the  Lagrangian  (A-22)  to  the  nonlinear 
hydrodynamics  can  be  made  by  replacing 


{**/«* 


dx 

ax 


Pe 

p 


x(x+a)  -  x(x)  +  a  ♦  •••  , 

ax 

where  x(x,t)  is  the  location  of  a  particle  at  time  t  which  was  at 
point  x  at  time  "0",  so  that 


"  ■  (ffU 
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APPENDIX  B 


THE  EFFECT  OF  DECAY  BANDWIDTH  ON  ATTENUATION 

In  che  second  chapter,  we  calculated  the  attenuation  of  sound  due 
to  nonlinear  Interaction  with  a  background  and  showed  that  the  magni¬ 
tude  of  the  attenuation  was  sensitive  to  the  nature  of  the  dispersion. 
In  this  appendix,  we  show  that  this  sensitivity  is  absent  when  either 
Che  sound  or  background  waves  are  attenuating  significantly.  We  will 
also  calculate  the  change  in  the  wave  speed  that  results  from  the  non¬ 
linear  interactions. 

A  sound  wave  that  is  attenuating  has  an  inherent  uncertainty  in 
frequency,  which  is  called  the  decay  bandwidth.  If  this  uncertainty, 
which  is  dissipative  dispersion,  is  much  larger  than  the  reversible 
dispersion  then  the  net  effect  of  dispersion  on  the  attenuation  is 
negligible. 

To  calculate  these  nonlinear  effects,  we  use  a  perturbation  method 

27 

developed  by  Bedeaux  and  Mazur.  The  method  is  to  add  to  the  fluid 
equations  an  external  driving  force  and  a  fluctuating  force.  As  an 
assumption,  both  forces  drive  the  equations  in  an  approximately  linear 
manner.  The  purpose  of  the  external  force  is  to  identify  a  Green's 
function,  while  the  fluctuating  force  by  design  creates  the  background. 
The  presence  of  the  background  changes  the  Green's  function  through 
nonlinear  interactions.  From  this  change,  we  directly  calculate  the 
attenuation  and  speed  of  a  sound  wave. 

The  effect  of  the  decay  bandwldths  will  be  simulated  by  adding 


56 


to  the  fluid  equations  a  dissipative  stress  tensor  of  a  Navier-Stokes 
form.  The  attenuation  of  sound  caused  directly  by  this  addition  will 
be  ignored  but  the  resulting  bandwidth  which  will  appear  in  the  non¬ 
linear  attenuation  will  be  replaced  by  the  true  bandwidth. 

Before  applying  the  perturbation  method,  we  make  some  simplifying 
modifications  and  assumptions.  First  since  the  external  and  fluctu¬ 
ating  forces  are  weak,  we  expand  the  fluid  equations  in  the  small 
deviations  about  the  rest  fluid,  and  keep  terms  which  are  at  most 
bilinear  in  the  deviations.  Second  we  neglect  the  nonlinear  terms  in 
the  fluid  velocity.  (In  Chapter  II  ,  the  velocity  contributions  were 
calculated  by  a  different  method.)  Third  we  assume  that  the  wavenumber 
of  the  sound  wave  is  much  less  than  the  background  wavenumbers. 

Finally,  we  will  assume  that  the  dissipative  and  reversible  dispersions 
are  small. 

In  Chapter  I,  we  derived  the  dispersive  fluid  equations: 
hs  ♦  T 

r-  ♦  vvs  ■  0  (conservation  of  entropy)  ,  (B-l) 


—  +  7»(pv)  »  0 
dt 


(conservation  of  mass) 


and 


(B-2) 


£ ,  22_1 

at  o  a,  l  ap  3x.a7.pJ 

To  expand  these  fluid  equations  in  small 
fluid,  we  let 


(Newton's  law)  .  (B-3) 

deviations  about  the  rest 


p  •  p  +  6p(x,t)  ,  v  •  6v(x,t)  ,  and  s  ■  s  ♦  5s(x,t)  , 

e  € 

where  and  s^  are  the  density  and  specific  entropy  of  the  rest  fluid 

and  the  quantities  with  the  prefix  6  are  the  deviations  from  rest.  We 
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then  substitute  these  expressions  into  the  fluid  aquations,  Taylor 
expand  the  internal  energy  about  the  rest  state,  and  keep  only  terms 
at  most  bilinear  in  the  deviations.  We  begin  by  considering  the 
entropy  equation  (B-l) : 

—  +  v*7s  ■  0 
3t 

Substituting  the  deviation  expressions,  we  find: 

&  -►  ** 

— — 6s  +  6v»75s  »  0 
at 

This  equation  states  that  the  deviation  in  entropy  is  proportional  to  a 
small  quantity  times  itself,  which  implies  that  the  entropy  deviation 
is  approximately  zero  and  entropy  per  mass  is  constant.  Since  entropy 
is  constant,  the  expansion  of  the  internal  energy  to  second  order  in 
the  density  deviation  is 


U(p,s,7p)  -  U(p  ,s  ,0)  ♦  SS. 


e'  e' 


*  ap2 


6p2  +  - ■’  (7p)2 


a($p)2 

and  the  two  remaining  fluid  equations  (B-2,3)  in  terms  of  the  devia¬ 
tions  are: 

— — +  7*(6p6v)  +  p  7»6v  ■  0  and 
3t  e 


— —  +  (6v*7)6v  - - 76p  -  —  7(6p)2  -  a7  726p  , 

dt  P  2p 


where 


a3ui 
ap3  ‘ 


e 


and 


au 


S3($p)2 


If  we  take  the  time  derivative  of  the  first  equation,  subtract  the 
gradient  of  the  second  equation,  and  neglect  nonlinear  terms  in  the 


S8 


velocity,  we  have  one  equation  for  mass  density: 

0  6p  -  -  -ic  726p2  ,  (B-4) 

where 

0  s  —  -  C2V2  +  cc74 
3t2 

If  we  had  included  the  Navier-Stokes  dissipative  tensor  in  the  fluid 
equation  pertaining  to  Newton's  law  (B— 3) ,  we  would  have  one  altera¬ 
tion: 

Q  i  -  c272  ♦  a74  ♦  v.  72  —  , 

8t2  b  3t 

where  is  four-thirds  the  shear  viscosity  plus  the  bulk  viscosity, 
which  will  be  assumed  to  be  a  constant. 

To  apply  the  Green's  function  method,  we  add  a  driving  force  and 
a  fluctuating  force  to  the  one  remaining  fluid  equation  (B-4) ,  which 
is  our  basic  equation  of  interest: 

0  4p  +  J  725p2  -  T(x, t)  +  F(x, t)  ( B— 5 ) 

where  T  is  the  fluctuating  force  and  F  is  the  driving  force.  Next  we 
Fourier  analyze  this  equation: 

(-  w2  *  W(k)2)  5 p  -  -  k26p2  -  T  ♦  F 
2 

2  2  2  2  a 

where  6p,  5p  ,  T,  and  F  are  now  Fourier  transforms  and  W(k)  He  k  +ak  + 

iv.GJk2. 

D 

When  the  forces  are  weak,  we  can  initially  neglect  the  nonlinear 
term,  and  treat  the  density  deviations  created  by  each  force  as  co¬ 
existing  independently.  In  this  approximation,  we  let  5o»6o_+6d_, 

r  T 

where  6p^,  and  SOj.  are  created  by  the  fluctuating  and  driving  forces, 
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respectively.  The  deviations,  dpF  and  SpT,  satisfy  the  equations: 

6pf  «  GqF  and  6pt  ■  GQT  ,  (B-6) 

-12  2 

where  G  is  the  linear  Green's  function:  G  » (— tu  +W(k)  ). 
o  o 

The  nonlinear  term  in  the  fluid  equation  (B-5)  is  the  source  of 

interaction  between  <Sp_  and  <5p_.  Including  this  term,  we  find  a  more 

T  r 

accurate  equation  for  5pF: 

(G"1  f  L)&PF  -  F 

where  L  is  the  following  operator: 

L  -  -  eic26pT  , 

where  it  is  implicit  that  L  involves  a  convolution: 

Lf  ■  -  e2k2[dq  6pT(k-q) • f(q)  . 

We  then  formally  invert  the  equation  for  Spy  and  ensemble  average  over 
the  fluctuating  quantity  6p^: 

6pF  -  GF  ,  (B-7) 

where  G,  which  is  identified  as  the  Green's  function  in  the  presence 
of  the  background,  is 

G  ■  <(L  ♦  G0L)"1>G0  . 

From  the  singularities  of  the  Green's  function,  we  will  determine 
both  the  attenuation  and  speed  of  a  sound  wave.  To  do  this  we  let  the 
external  driving  force  be 


F(x, t) 


exp( ik*r)6(t) 
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and  chen  use  Che  Green's  funcCion  and  Eq.  (B-7)  to  find  Che  response 


<50 


F* 


fiPi 


exp< 


ik*r)|  i 


dw  G(w,k)  exp(iwt) 


Using  contour  Integration,  we  Identify  the  propagating  nodes  for  wave- 
number,  k,  to  be 


6p  -  exp(ik*r  ♦  iu.t) 
F  J 


where  are  the  simple  poles  of  the  Green's  function  for  lmu>0.  The 
The  Imaginary  part  of  is  the  attenuation  coefficient  and  the  real 
part  divided  by  the  wavenumber  is  the  wave  speed. 

For  our  fluid  equation,  the  Green's  function  will  be  of  the  form: 

6  -  [  u2  -  W(k)2  ♦  g(w,k)  ]“l  , 

where  che  function  g  will  be  transcendental.  To  find  the  singularities 
of  the  Green's  function,  we  will  presume  that  the  presence  of  function 
g  does  not  significantly  alter  the  positions  of  the  poles  so  that  the 
dependence  of  g  on  u  can  be  removed  by  setting  u>  equal  to  its  approxi¬ 
mate  value  in  the  vicinity  of  the  pole: 

G  -  [  u2  -  W(k)2  ♦  g(W(k),k)  ]"1  , 

The  pole  of  the  Green's  function  is  at  the  location: 

_  l_ 

uj  -  [  W(lt)  -  g(W(k),k)  ]2 

Expanding  the  square  root  in  a  Taylor  series  about  W(k),  we  find: 


wj  "  -  -j  g(W(k) ,k)/W(k) 


61 


The  imaginary  pare  of  uj^  is  Che  attenuation  coefficient,  B: 


3  -  la  g(W(k),k)/2W(k)  , 


(B-8) 


and  the  real  part  divided  by  the  wavenumber  is  the  wave  speed  c (k) : 

c(k)  -  Re  M  W(k)  -  i  g(W(k)  ,k)/W(k)  ]  (B-9) 

k  *  • 

To  find  the  Green's  function  in  the  presence  of  the  background 

(6ot),  we  must  evaluate  the  ensemble  average  of  the  quantity:  (l+GoL)-1 

We  average  this  quantity  by  first  expanding  in  powers  of  L,  which  is 

presumed  small,  and  then  approximating  by  keeping  the  first  three 

terms.  The  terms  are  1,  GqL,  and  GqLGqL.  The  average  of  the  first 

term  is  itself.  For  the  second,  if  we  impose  translational  invariance, 

then  <5p  >«0  and  since  -G  L-G  ek  <5p_,  we  deduce  that  <-G  L>«0.  For  the 
1  O  O  T  o 

third  term,  we  write  the  term  acting  on  the  driving  force  with  explicit 
convolutions: 

GQLG0LF  -  e2k2  Gq  (<*>,k)  JJJ J  d«'d3k'd«"d3k"C0  (w",k")k,,2F(u>'  ,k' )  * 


6p^(k"-k'  , k” 

The  ensemble  dependence  in  this  expression  is  in  the  quantity: 

*  *  *  * 

6pT(k"-k' ,u"-<i),)-6pT(k-  k'\u>-<D") 

whose  ensemble  average  requires  knowledge  of  the  fluctuating  force. 

To  proceed  we  must  design  the  fluctuating  force  to  create  our 
desired  background  which  we  assume  is  isotropic  with  spectral 
intensity  u(k).  To  meet  this  demand,  we  let  the  fluctuating  force 
satisfy  the  following  correlation  relation: 
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<T(X,q)T(X*  ,q'  )>  -  -e  .b-  q2u(q)  6(X-*-X')  6(q+q')  . 

4* 

One  can  Chen  verify  that  for  this  correlation  of  Che  fluctuating  forces 
Che  average  energy  density  per  unit  wavenumber  of  the  sound  field  is 
given  by 

[jj  dk'dw'du  <6pT(a},ic)6pT(u>,  )>  »  u(k)  , 

where  (B-6)  is  used  for  6p^,.  Using  this  correlation  relation  and  the 
equation  relating  6pT  and  T,  we  have  the  necessary  ensemble  average: 

<6pt(x,q)6pT(x',q')>  -  q2u(q)|G0(X,q)j26(x+X')  6(q+q'  )  . 

The  presence  of  the  delta  functions  is  a  consequence  of  imposing  trans¬ 
lational  Invariance  for  space  and  time. 

We  have  shown  that  the  Green's  function  in  the  presence  of  a 
background  is 

G  »  (  1  ♦  G0g(u>,10  )GQ 

where 

g(w,ic)  =  <LGqL>  . 

Since  the  function  g  is  small,  we  rewrite  this  equation  in  the  approxi¬ 
mate  form: 

G  ■  (  G'1  ♦  gCu.t)  . 

When  we  neglect  the  viscosity  term  in  G  ,  we  have  shown  that  the  atten- 

o 

uation  coefficient  (S)  and  the  wave  speed  are  simply  related  to  the 
function  g  by  Eqs.  (B-8,9). 
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What  remains  to  be  done  is  the  evaluation  of  the  function 


gF  -  e2k2J[[|  du),d3k'dw"d3k,,GQ  (w" ^"Jk"2?^1  ,k‘ )  * 

— m 

<6pT(k,,-k\<tf,-w,)*6p  (k-  k",W(k )-<*>")> 

Substituting  the  expression  for  the  ensemble  average  of  squared  density 
deviation,  making  the  change  of  variables  K»k-k"  and  r*W(k)-w",  and 
integrating  the  single  primed  integrations,  we  find: 

8  *  “  e2Pe  k2fd3K  u(R)(K-k)2K2[dr  G0(r,K)G0(-r,-K)G0(W(k)-r,k-R)  . 

Next  we  use  contour  integration  to  perform  the  integral  over  T,  which 
gives  the  result: 

8  "  ”  ~2E2Pe^fd3Ku(KKK-k)2  I[(W(k)±W(K)+^(K2+(k-K)2))2-W(k»K)^'1  • 

to  integrate  over  R  space,  we  will  use  spherical  coordinates  whose 
x-axis  is  chosen  along  k.  The  integration  over  azimuthal  angle  is 
trivial  and  has  Che  value  twice  pi. 

Before  proceeding  to  perform  the  integration  over  the  polar  angle, 
we  will  approximate  the  denominators  of  the  integrand.  Since  the  back¬ 
ground  wavenumbers  are  assumed  to  be  much  greater  than  the  wavenumber 
of  the  sound  wave,  we  will  keep  only  the  leading  order  in  the  ratio 
k/R.  The  dispersion  is  assumed  small  so  that  we  keep  only  the  leading 
order  in  a  and  With  these  approximations,  the  expression  for  g 

becomes 


64 


65 


Before  finding  the  attenuation  coefficient  ana  -ive  speed,  we  make 


three  revisions.  First  we  remove  the  dependence  of  g  on  viscosity  by 

2 

replacing  the  quantity  - =•  ,  which  is.  the  lifetime  of  a  wave  of  wave- 

v  KZ 

number  K  in  the  viscous  °fluid,  by  the  true  lifetime  which  we  will 

label  t(K) .  Second,  we  include  the  effect  of  nonlinear  terms  in 

velocity  by  redefining  £  as  follows: 

,  s  2c2. 
e--A 


where  as  before 

.  ,  p  dc 

c  dp 

Finally  we  assume  that  the  logarithm  and  arctangent  in  the  expression 
for  g  do  not  vary  much  over  the  range  of  wavenumbers  for  which  the 
background  energy  density  (u(k))  is  significant.  This  assumption 
allows  us  to  take  the  logarithm  and  arctangent  outside  the  integral  and 
evaluate  them  at  the  wavenumber  for  which  the  background  enery  is 
peaked.  Incorporating  these  revisions,  we  find: 


I  A2k2E  J  2  +  ylnf  - —  - 

Pe  [  1  +  ((3aK2/4c2)aisT)2 

i|  arctan(u)9t)  +  arctan((3aK2/4c2)u3x) 


where  K  is  the  wavenumber  for  u(k)  is  maximum,  E  is  the  total  energy 
density  (i.e.,  E"/dku(K)),  and  t»t(K). 

The  relationships  between  the  function  g  and  the  attenuation  coef¬ 
ficient  and  wave  speed  was  shown  to  be  Eq.  (B-8,9).  Using  these  rela- 

28 

tionships  and  the  calculated  g,  we  find  to  leading  order  in  a  and  v  : 

D 


p  - 


A2 

pec' 


“.E( 


arctan{ojst)  +  arctan((3aK 2/4c2)usr) 


and 


c(u>s) 


c0  + 


A2 
2PeC0 


E<ln 


1  +  (u>sr)2 


1  ♦  ((3<xK2/4c2)uist)2 


+  4 


where  cq  Is  the  wave  speed  in  the  absence  of  nonlinearities  or  thermal 
mo cion. 

The  assumptions  and  approximations  that  we  have  made  in  calculat¬ 
ing  Che  wave  speed  and  attenuation  are  appropriate  to  describe  sound 
waves  in  liquid  helium  at  low  temperatures.  The  background  is  thermal 
motion  for  which  we  know  from  statistical  mechanics  the  total  energy 
density: 

4«5  (kbT)*» 

E  «  ■- 

15h3c3 

where  T  is  temperature,  kg  is  Boltzmann's  constant,  and  h  is  Planck's 
constant. 

The  attenuation  and  wave  speed  depend  on  the  product  of  the  sound 

frequency  and  the  lifetime  of  the  background  waves.  When  this  product 

is  much  smaller  than  one,  the  attenuation  coefficient  becomes  propor- 

A2E  2 

tional  to  the  sound  frequency  squared:  6*  - 0  tui  ,  and  the  change 

Pc-  s 
0 

in  wave  speed  due  to  nonlinear  effects  is  frequency  Independent: 

2A2E  2 

.  When  the  product  (<u  t)  lies  in  the  range  1<<w  r<<  — ■  -  , 

e  8  3oK2 

then  the  attenuation  coefficient  is  the  same  as  Westervelt's  result  for 

a  nondispersive  fluid,  and  the  change  in  wave  speed  is 
A2  7 

Ac  *  7TZ  E  (4  +  ))  • 
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When  the  product  oj  t  is  large  so  that  — —  u)  T»l,  then  the  attenuation 
s  z  s 

c 

coefficient  is  either  zero  or  twice  Westervelt's  result  depending  on 


whether  the  dispersive  parameter  a  is  negative  or  positive,  respectively. 

2ajr2 

The  change  in  wave  speed  for  u  T  — r—  >>1  becomes  again  frequency  inde- 

a2_  t  2  >2}  8  c 

pendent:  Ac  *vrr-(4+ln  “ — r  •  For  liquid  helium  at  0.0  bars  and  .3® 

pe  '•3aK,iJ  2 

Kelvin,  the  value  of  the  quantity  — c--_  is  predicted  by  experiment5  to 

3  ax 

be  roughly  one  hundred  if  the  dispersion  relation  is  assumed  to  be  of 
the  form  (1-34) . 


In  liquid  helium,  neither  the  lifetime  of  the  background  waves  nor 
the  fine  details  of  the  dispersion  have  been  determined.  For  our  re¬ 
sults,  it  is  clear  that  no  value  of  the  lifetime  completely  accounts 
for  the  experimental  results.  The  discrepancy  may  be  due  to  our 
implicit  assumption  that  the  sound  wave  lifetime  was  longer  than  that 
of  the  background  waves. 
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APPENDIX  C 


WIEN'S  DISPLACEMENT  LAW  FOR  A  PONDERABLE  CONTINUUM 

29 

In  the  late  nineteenth  century,  Wien  derived  the  displacement 
law  which  states  that  the  electromagnetic  radiation  in  equilibrium  has 
its  energy  distributed  over  wavelengths  according  to  the  form: 

u(k)  -  k3  g(k/T)  , 

where  u  is  the  energy  per  unit  wavenumber,  T  is  temperature,  k  is  wave- 
number,  and  g  is  an  undetermined  function.  For  a  fluid,  we  will  prove 
that  the  thermal  energy  is  distributed  over  wavelengths  in  the  follow¬ 
ing  way: 

u(k)  •  ck3  f(ck/T)  .  (C-0) 

where  c  is  the  speed  of  sound  and  f  is  an  undetermined  function.  The 

appearance  of  the  speed  of  sound  is  due  to  the  fact  that  the  sound 

speed  is  not  constant,  unlike  the  speed  of  light. 

Our  derivation  of  the  displacement  law  for  a  fluid  is  outlined  as 

30 

follows.  First  we  show  following  the  method  of  Whitham  that  for  a 
sound  wave  the  wave  energy  divided  by  the  frequency  is  an  adiabatic 
invariant.  Next  we  consider  a  quasistatic  change  in  volume  and  use  the 
fact  that  the  ratio  of  wave  energy  to  frequency  is  constant  for  each 
thermal  component  to  calculate  the  change  in  thermal  distribution. 

From  this  change,  we  will  show  that  the  thermal  distribution  as  a 
function  of  the  volume  and  wavenumber  has  the  form: 


u(k)  -  ck3  h(k3V) 


where  V  is  Che  volume  and  h  is  an  undetermined  function.  We  then  use 

3 

VT 

thermodynamics  and  hydrodynamics  to  show  that  the  quantity  — remains 

c 

constant  for  a  quasistatic  change  in  volume.  Using  this  fact  and  the 
functional  dependence  of  the  thermal  distribution  on  volume,  we  deduce 
dependence  of  the  thermal  distribution  on  temperature  is  of  the  form 
(C-0) . 

That  the  ratio  of  the  wave  energy  to  frequency  Is  an  adiabatic 
Invariant  will  be  shown  using  the  principle  of  least  action.  Action  is 
defined  as  the  Integral  of  the  Lagrangian  density  which  is  the  dif¬ 
ference  between  the  kinetic  and  potential  energy  densities.  To  de¬ 
scribe  a  sound  wave,  we  will  show  that  the  Lagrangian  can  be  approxi¬ 
mated  and  written  in  terms  of  the  independent  variables  relating  to 
the  amplitude  and  phase  of  the  wave.  Demanding  that  action  be  a 
minimum,  we  will  show  that  resulting  dynamical  equations  imply  that  the 
ratio  of  energy  to  frequency  is  an  adiabatic  invariant. 

The  Lagrangian  density  that  is  appropriate  for  a  sound  wave  in  a 
fluid  is  of  the  form: 

L  »  pv2  -  U(p)  , 

where  p  is  mass  density,  v  is  the  fluid  velocity,  and  U  is  the  internal 
energy.  Since  the  sound  wave  is  presumed  to  be  of  low  amplitude,  we 
expand  the  Lagrangian  density  in  terms  of  the  deviations  and  keep  only 
leading  order  terms: 

L  .  i  p^*2  .  U(p  )  -  S£|  5p  -  6p2  , 

op  I e  2p 
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where  Pe  is  Che  rest  density,  <$v  and  6p  are  the  deviations,  and 
2  ^2  U 

c =  •  The  leading  term  in  the  expansion  of  the  internal  energy 

is  neglected  because  it  is  constant  and  will  have  no  effect  when  de¬ 
manding  that  action  be  minimum.  The  second  term  is  also  neglected  be¬ 
cause  its  contribution  to  the  sound  wave  energy  is  zero.  The  effective 
Lagrangian  density  is  therefore: 

L  ■  i  p  Sv2  -  —  6p2 
z  e  2p 

e 

Next  we  rewrite  the  Lagrangian  in  terms  of  the  Lagrange  coordinates 
using  the  relations: 

♦  •*  . 

v  "  at  Md  6p  *  pe7*n 

where  n  is  the  displacement  of  the  fluid  particle  from  equilibrium. 

For  a  sound  wave,  the  displacement  can  be  written  in  terms  of  the 
amplitude  and  phase: 

♦  + 

h  -  A  exp(i9)  +  c.c.  , 


where  A  is  the  amplitude  and  9  is  the  phase  to  which  frequency  and 

wavenumber  are  related  as  follows: 

as  *  * 

ij>  m - and  k  *  7  9 

at 

When  the  density  of  the  fluid  is  changing  slowly,  then  the  quantities 
A,  to,  and  k  are  also  slowly  varying.  Substituting  the  form  of  the 
displacement  into  the  Lagrangian  density  and  neglecting  the  slow  vari¬ 
ations,  we  find 

L  •  -  Re[AexpU9)  -  Aexp(-i9)  ]2{  (||) *  ♦  c2(79)2^ 


Since  we  are  interested  in  global  properties  of  the  sound  wave,  we 
average  the  Lagrangian  over  a  wavelength: 


<L>  “  Pel^|2  (t||)2  +  =2^®)2) 

Performing  the  variations  of  the  independent  variables  A  and  6,  and 
demanding  that  action  be  a  minimum,  we  find  the  equations: 

Iff)  *  =2<70)2  “d  i7  (Pe!A!2i7  ]  *  7(PelA|2c2v  e)  ■  0  • 

Using  vector  analysis  and  algebra,  we  can  use  the  two  equations  above 
to  derive  the  following  equation: 

h q  *  H  ifr cQ  1  ■ 0  “here  Q  ■  p«i*i2{(I!1!*  c2<;-9)2}  *  (0 

This  equation  is  in  the  form  of  a  conservation  equation  for  which  the 
invariant  quantity  is  Q.  The  excess  energy  density  for  the  fluid  is 


approximately: 


E  ■  p  6v^  *  -2—  <5p^ 


which  if  put  in  terms  of  A  and  9,  and  averaged  gives  the  energy  of  the 


sound  wave: 


v  a© 
‘'aT”  • 


Since  the  quantity  yj-  is  the  frequency  and  the  quantity  Q  is  invariant, 
we  deduce  that  the  energy  divided  by  frequency  is  an  adiabatic  invari¬ 
ant.  For  the  particular  case  of  the  quasistatic  change  in  volume  we 


—  (E/w)  »  0 
dV 


(C-l ) 


We  now  consider  a  small  quasistatic  change  in  volume  of  a  fluid 
at  some  initial  temperature.  The  energy  of  thermal  motion  can  be 
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viewed  as  being  distributed  in  an  unique  way  over  the  wavenumbers  of 
che  thermal  waves.  The  quasistatic  change  in  volume  alters  the  distri¬ 
bution.  If  nonlinear  interaction  between  waves  is  neglec cable,  then 
the  change  in  the  distribution  can  be  determined  by  applying  (C-l)  to 
each  of  the  waves  chat  compose  the  thermal  distribution.  From  the 
change  in  distribution  we  will  derive  information  about  the  form  of  the 
distribution. 

The  proof  that  nonlinearities  are  negiectable  is  as  follows.  We 
imagine  that  the  fluid  is  compressed  with  nonlinearities  turned  off 
so  that  the  distribution  changes  according  to  (C-l) .  Now  we  turn  on 
the  nonlinearities.  If  we  started  from  an  equilibrium  state  and  if  the 
compression  is  reversible  then  there  can  be  no  entropy  change  due  to 
nonlinearities  and  the  new  distribution  generated  from  (C-l)  must  be  an 
equilibrium  distribution  for  which  nonlinearities  have  no  effect. 

To  find  the  change  in  distribution  due  to  the  quasistatic  volume 
change,  we  let  h  be  the  mechanical  function  [h«h(k,7)]  such  that  the 
energy  density  in  the  wavenumber  interval  dk  at  k  is  hck  dk  (i.e.  we 
start  with  a  certain  distribution  h  over  k  for  the  given  volume  and 
now  change  volume  and  therefore  take  h»h(k,V)).  Then  the  fact  that  the 
energy  in  dk  divided  by  its  frequency  (ck)  is  invariant  implies  the 
relation: 

Vhk2  dk  -  V'h'k'2  dk'  ,  (C_2) 

where  V  is  the  volume  and  the  primes  denote  the  changed  quantities. 

This  relation  can  be  simplified  by  expressing  the  primed  wavenumber  in 
terms  of  the  change  in  volume.  Boundary  conditions  at  the  walls  con¬ 
taining  the  fluid  demand  that  the  wa/enumber  be  inversely  proportional 
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to  the  linear  dimension  of  the  system  and  since  volume  is  proportional  to 
Che  third  power  of  the  linear  dimension,  we  can  use  simple  calculus  to 
derive  the  following  relation  which  is  valid  to  first  order  in  the 
change  in  volume: 

k'  "  (  1  "  37  1  k  *  (c"3) 

where  AV  is  the  change  in  volume.  Substituting  the  k'  from  (03)  into 
(02)  and  keeping  terms  to  order  fit,  we  find  h*h\  (C-4) 

Using  Eq.  (04),  we  will  show  that  the  thermal  distribution  is  of 
the  form: 

u(k)  *  ck3h(k3V)  .  (05) 

First  we  expand  the  primed  h  in  a  Taylor  series  about  the  initial 
wavenumber  and  volume,  and  keep  only  terms  first  order  in  AV: 

h'  -  h(k\V')  -  h(k, V)  +  Ak  +  |§AV  .  (06) 

Using  the  Eqs.  (03),  (04),  and  (06),  we  can  relate  the  partial 
derivatives  of  h:  y  •  V  |^  .  The  solution  to  this  partial  differ¬ 
ential  equation  is  h«h(k3V),  which  proves  the  thermal  distribution  is 
of  the  form  (05). 

Our  ultimate  goal  is  to  derive  the  displacement  law,  which  will 

3  3 

follow  from  (05)  and  the  fact  chat  the  quantity  VT  /c  is  a  tfermo- 
dynamic  adiabatic  invariant.  The  proof  that  VT3  /c3  is  invariant  will 
proceed  as  follows.  First  we  separate  both  the  pressure  and  energy  of 
the  fluid  into  two  pieces:  one  piece  due  to  thermal  motion  and  Che 
other  being  that  which  exists  in  the  absence  of  thermal  motion  (at 
zero  degrees).  Then  we  use  (C-l)  to  derive  the  relationship  between 
the  pressure  and  energy  of  the  excitations.  With  this  relationship 
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A  simple  relationship  between  and  can  be  derived  from  (C-l)  , 
(C-10),  and  the  assumption  that  the  speed  of  sound  depends  only  on 
density.  For  a  single  wave,  (C-l)  implies  that  in  a  quasistatic  change 
in  volume  the  change  in  wave  energy  in  terms  of  the  change  in  frequency 
is  given  by  the  equation: 

AE  _  Au) 

E  cj 

Since  oj-ck,  the  change  in  energy  can  be  expressed  in  terms  of  Ak  and 
Ac : 
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(C-ll) 


Ag  m  Ac_  +  Ak 


If  c  i3  only  a  function  of  density,  then  the  change  in  c  in  terms  of 


AV  is: 


ac  av 

8p  V 


(C-12) 


Using  (C-3),  (C-ll)  and  (C-12)  we  can  express  AE  in  terms  of  AV: 


(3  c  dp  >  V 


Since  this  equation  is  true  for  each  wave  that  composes  the  thermal 
distribution,  it  is  also  true  for  the  whole  distribution  so  that 


MUV)  ■  -  (  5  *  f  %  • 


Comparing  this  equation  to  (C-10) ,  we  deduce  that  P  is 

x 


p  -  +  u 

x  V  3  c  dP  J  x 


(C-13) 


As  the  relation  for  AE  also  applied  only  to  adiabatic  transforma¬ 


tions,  the  dependence  of  on  T  can  be  found  using  (C-13)  and  the 

second  law  of  thermodynamics.  A. consequence  of  the  second  law,  which 

arises  from  equating  mixed  partial  derivatives  of  entropy,  is  the 

. _.31 


equation: 


3UV  Ap 

av  t  ax  v  P 


(C-14) 


Since 


P„  ,  we  can  write  (C-1A)  for  just  P  and  U  : 
o  XX 


av  t 


3p, 

T  —  -  ? 
3T  x 


Substituting  the  P  from  (C-13),  we  have: 
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( -  * 

P  dc  | 

r .  a°xi  . ..  i 

H 

► 

<o 

l  3 

c  dp  J 

si 

< 

c 

X 

If  we  lee  Z“V/c 3,  we  can  transform  Che  last  equation: 

a _ to,  _  a _  vuy 

ainz  T  ainT3  T 

The  general  solution  to  this  equation  is 

0  -  -  FdnCzT3) 
x  7 

where  F  is  as  vet  undetermined.  Since  U  is  intensive,  the  function  F 

x 

must  be  exponential  and  therefore  the  total  energy  density  is  given  by 


constant  T4/c3 


CC- 15) 


Using  (C-10)  and  (C-15) ,  we  can  prove  that  the  quantity  VT3/c3  is 
invariant  for  a  quasistatic  change  in  volume.  Substituting  the  from 
(C-13)  and  the  from  (C-15)  into  (C-10),  we  find  with  some  algebra: 

aTl  .  f  1  *  P  dc  I  T  , 

avis  1  ^  c  dP  I  V  ’  (C-16) 

The  identity 


_1_  +  £_  dc_  m  c3  5  V 
3  c  dp  3  avc3 


is  useful  for  the  derivation, 
directly  that 


a_ 

av 


Using  equation  (C-16) ,  we  can  show 


-  0 


(C-17) 


An  implication  that  can  be  made  from  (C-J7)  is  that  during  any 

quasistatic  change  in  volume,  the  volume  is  proportional  to  a  constant 
3  3 

times  T  /c  .  Substituting  this  form  of  V  into  (C-5) ,  we  have  the 
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displacOMac  law  for  a  fluid: 


u(k)  ■  ck3  f (ck/T) 


rs 
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